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GENE DISCOVERY THROUGH COMPARISONS OF NETWORKS 
OF STRUCTURAL AND FUNCTIONAL RELATIONSHIPS 
AMONG KNOWN GENES AND PROTEINS 

SPECIFICATION 
The invention described herein was funded in part by a grant from the 
National Library of Medicine, namely, Grant Number's LM06274 and LM05627. 
The United States Government may have certain rights to the invention. The present 
specification contains a computer program listing which appears as a microfiche 
Appendix H. 

STATEMENT REGARDING MATERIAL SUBJECT TO COPYRIGHT 
A portion of the disclosure of this patent document contains material 
which is subject to copyright protection. The copyright owner has no objection to the 
facsimile reproduction by anyone of any portion of the patent document, as it appears 
in any patent granted from the present application or in the Patent and Trademark 
Office file or records available to the public, but otherwise reserves all copyright 
rights whatsoever. 

An appendix containing source code listing utilized in practicing an 
exemplary embodiment of the invention is included as part of the Specification, 



1. INTRODUCTION 

20 The present invention relates to methods for identifying novel genes 

comprising: (i) generating one or more specialized databases containing information 
on gene/protein structure, function and/or regulatory interactions; and (ii) searching 
the specialized databases for homology or for a particular motif and thereby 
identifying a putative novel gene of interest. The invention may further comprise 

25 performing simulation and hypothesis testing to identify or confirm that the putative 
gene is a novel gene of interest. 

The present invention relates to natural language processing and 
extraction of relational information associated with genes and proteins that are found 
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in genomics journal articles. To enable access to information in textual form, the 
natural language processing system of the present invention provides a method for 
extracting and structuring information found in the literature in a form appropriate for 
subsequent applications. Specifically, the present invention provides for the 
5 generation of specialized databases containing information on gene/protein structure, 
function and regulatory interactions based on the retrieval of such information from 
research articles and databases, and computer representation of such information in a 
manner that allows efficient access to the extracted information. 

The invention further provides for the use of the specialized databases 

10 for identifying novel genes based on detection of sequence similarities and 

domain/motif matches between genes/proteins, computation and interpretation of 
phyiogenetic trees for multigene families, and analysis of homologous regulatory 
networks. The methods of the invention are based on the observation that functionally 
similar regulatory systems are generated during evolution by genetic duplication of 

1 5 ancestral genes. Thus, a comparison of homologous/similar networks within the same 
organism and between different species will allow the identification of genes absent in 
one of the systems under comparison. In this way genes that contribute to the 
phenotype of a specific disease associated with a particular biological system under 
analysis may be identified. 

20 2. BACKGROUND OF THE INVENTION 

2.1. NATURAL LANGUAGE PROCESSING 
Researchers working in molecular biology must constantly consider the 
information present in the literature relating to their regulatory systems of interest and 
the genes and proteins that operate within those systems. Unfortunately, to remain up- 

25 to-date on the relevant literature, the researcher is required to perform laborious 
reading and manual integration of research articles, each of which may address a 
narrow subject. Therefore, technology that enables rapid retrieval of information from 
literature and manipulation of derived functional data should have a dramatic effect on 
the accesss of the researcher to important facts and ultimately should facilitate the 

30 discovery of novel human genes. 
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Natural language processing is an automated system that provides for 
a complex of programs for automatic retrieval of information from text analysis and 
for the computer representation of that information in a form that allows efficient 
access and extraction of that information. MedLee (Medical Language Extraction and 
5 Encoding System) has recently been successfully used for processing different types 
of medical texts as described in co-pending United States Patent Application Serial 
Number 09/370,329, incorporated herein in its entirety by reference (see also, 
Friedman et al., 1994, J. Amer. Med. Inf. Assoc. 1:161-174; Hripcsak et al. 1995, 
Ann. Intern. Med. 122:681-688; Hripcsak et al., 1998, Meth. Inform. Med.; Jain et al., 

10 1996, Proc. AMIA Annu. Fall Symp. 542-546; Knirsch et al, 1998). When tested, 
MedLEE was on average as successful in retrieving reports associated with specified 
clinical connections as twelve medical experts invited for evaluation of the system. 

Another text analysis technique has recently been developed that 
combines finite-state machines with statistical machine learning approaches. These 

1 5 models extract detailed semantic information from texts (e.g., see Hatzivassiloglou 
1996, In Klavens, J.L., and Resnick, P.S. (eds) The Balancing Act: Combining 
Symbolic and Statistical Approaches to Language, MIT Press, Cambridge, MA) when 
extensive prior knowledge about the domain is not available. The techniques have 
been subsequently applied to the tasks of (i) automatically identifying medical terms 

20 for the automated summarization of research articles reporting on clinical studies and 
(ii) sanitizing sensitive information in patient records so that they can be widely 
disseminated for research purposes. 

A number of projects have also been developed as statistical 
information extraction tools that operate with limited or no prior knowledge about the 

25 application domain. These earlier efforts include XTRACT, a tool that recovers 
collocational restrictions between words that has been licensed to more than thirty 
sites worldwide (Smadja, F., 1993, J. Comp. Ling. 19:143-177), CHAMPOLLION, a 
system that retrieves bilingual mappings between words and phrases in parallel texts 
from different languages (Smadja, F. et al. 1996, J. Computational Linguistics 22:1- 

30 38), and a system that automatically aligns noisy, semi-parallel texts from different 
languages (Fung, P. and McKeown, K.R., 1997, Machine Translation 1 1 :23-29). 
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2.2. IDENTIFICATION OF NOVEL GENES 
A variety of different methods are currently utilized for the 
identification and characterization of novel genes. Perhaps the most widely used 
5 method for generating large quantities of sequence information is via high throughput 
nucleotide sequencing of random DNA fragments. A disadvantage associated with 
this gene discovery technique is that in most instances when genes are identified their 
function is unknown. 

For identification of specific disease genes, positional cloning is 

10 currently the most widely used method. The positional cloning approach combines 

methods of formal genetics, physical mapping and mutation analysis and usually starts 
with a precise description of the disease phenotype and a tracing of the disease 
through families of affected individuals. Genetic linkage data obtained from the 
analysis of affected families frequently allows the determination of an approximate 

1 5 genomic localization of the candidate disease gene with a precision of several millions 
of nucleotides. Once localized, the genetically defined chromosomal region is then 
recovered from genomic libraries as a contiguous set of genomic fragments. Genes 
residing in the disease-related region are determined by analysis of transcripts that are 
transcribed from the genomic fragment. From this analysis an initial set of candidate 

20 genes for a particular disease are identified based on the presence of the gene product 
in the biological system affected by disease and a correlation between its expression 
pattern and the pattern of disease progression. 

Important information for selection of candidate genes also comes from 
analysis of their homology with genes known to be part of the same or related 

25 biological system. Finally, the ultimate proof of association between a gene and a 

genetic disorder comes from mutational analysis of a gene in patients affected by the 
disorder and from demonstration of a statistical correlation between occurrence of 
mutation and the disease phenotype. 

Although positional cloning is a powerful method for gene discovery, 

30 the experimental method is extremely tedious and expensive. Moreover, disease 

genes implicated in genetically complex disorders, i.e., those controlled by multiple 
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loci, can hardly be found using this strategy because of the complications associated 
with multiple loci linkage analysis. 

Specialized databases for homology searches have also been utilized in 
disease gene discovery projects. In recent years a number of efficient sequence 
5 comparison tools have been developed such as the BLAST (Basic Local Alignment 
Search Tool) family of programs designed for comparison of a single "search 
sequence" with a database (see Altschul et al., 1990, J. MoL Biol. 215:403-410; 
Altschul et al., 1997, Nucleic Acids Res. 25:3389-3402), the family of Hidden 
Markov Model methods for comparison of a set of aligned sequences that usually 
10 represent a protein motif or domain with a database (e.g., Krogh et al., 1994, J. Mol. 
Biol. 235:1501-1531; Grundy et al., 1997, Biochem Biophys. Res. Commun. 231:760- 
6) and various other comparison tools (Wu et aL, 1996, Comput. Appl. Biosci 12:109- 
118; Neuwald et al., 1995, Protein Sci. 4:1618-1632; Neuwald, 1997, Nucleic Acids 
Res. 25:1665-1677). 

15 When used in disease gene discovery projects, homology searches can 

be enhanced by creating specialized databases that utilize statistical analysis for 
evaluating significance of sequence similarities in comparison of new sequences with 
a database of known sequence. Such databases are fine-tuned to the size of the 
database used (Altschul et al., 1990, J. Mol. Biol, 215:403-410; Altschul et al., 1997, 

20 Nucleic Acids Res. 25:3389-3402), so that the same level of homology between a 

search sequence and a database sequence can be determined to be highly significant if 
the search sequence is compared with a smaller database, or insignificant and thus 
undetectable, if the search sequence is compared with a larger database. 

In alternatives to standard homology searches, in projects oriented 

25 towards gene discovery, researchers usually have some a priori knowledge about the 
set of genes/proteins that might display important similarity to the unknown new gene. 
Therefore, selecting an a priori defined set of genes/proteins for comparison with new 
experimental sequences is a feasible and useful strategy. This strategy was 
successfully applied to search for homologs of disease genes in yeast and nematode 

30 genomes by Mushegian et al. (1997, Proc. Natl. Acad. Sci USA 94:5831-5836). 
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Two homologous genes taken from different species that originate 
from the nearest common ancestor by speciation are referred to as orthologs, while 
any two genes that originate from a common ancestor via a series of events involving 
intragenomic duplications are call paralogs. Tatusov et al. (1994, Proc, Natl, Acad. 
5 Sci USA 91 : 12091-12095) describe comparisons of proteins encoded by the genomes 
of different phylogenetic lineages and elucidation of consistent patterns of sequence 
similarities permitting the delineation of clusters of orthologous groups (COGs). Each 
COG consists of individual orthologous genes or orthologous groups of paralogs from 
different phylogenetic lineages* Since orthologs typically have the same function, the 

10 classification of known genes and proteins into clusters of orthologous groups permits 
the assignment of a function to a newly discovered gene or protein by merely 
classifying it into a COG. Although Tatusov describes a method for assigning a 
function to a newly discovered gene, he does not describe a method for predicting the 
existence of undiscovered genes. In addition, Yuan, et al. attempted simultaneous 

1 5 reconstruction of a species tree and identification of paralogous groups of sequences 
and detection of orthologs in sequence databases (Yuan et al, 1998, Bioinformatics 
143:285-289). 

Other groups have aimed at capturing interactions among molecules 
through the use of programs designed to compare structures and functions of proteins 

20 (Kazic 1994, In: Molecular Modeling: From Virtual Tools to Real Problems ^ 

Kumosinski, T. and Liebman, M.N. (Eds.), American Chemical Society, Washington, 
D.C pp. 486-494; Kazic, 1994, In: New Data Challenges in Our Information Age 
Glaesar, P.S. and Millward, M.T.L. (Eds.). Proceedings of the Thirteenth 
International CODATA Secretariat, Paris pp. C133-C140; Goto et al., 1997, Pac. 

25 Symp. Biocomput. p. 175-186; Bono et al, 1998, Genome Res. 8:203-210; Selkov 
et al., 1996, Nucleic Acids Res. 24:26-28). These projects are significantly different 
from the inventive methods described herein because they do not describe methods for 
deducing the existence of as yet unknown genes based on comparisons of regulatory 
pathways and gene structure between one or more species. The present invention 

30 provides a method for increasing the sensitivity of analysis methods through the 
generation of specialized databases. 
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3. SUMMARY OF THE INVENTION 
In accordance with the present invention there is provided methods for 
identification of novel genes comprising (i) generating one or more specialized 
databases containing information on gene/protein structure, function and/or regulatory 
5 interactions; and (ii) searching the specialized databases for homology or for a 
particular motif and thereby identifying a putative novel gene of interest. The 
invention may further comprise performing simulation and hypothesis testing to 
identify or confirm that the putative gene is a novel gene of interest. 

The invention is based, in part, on the observation that functionally 
10 similar regulatory systems are generated during evolution by genetic duplication of 
ancestral genes. Thus, by comparing phylogenetic trees or regulatory networks and 
identifying genes and/or proteins absent in one system under comparison, the 
existence of as yet unidentified genes and/or proteins can be predicted. To make 
meaningful comparisons of phylogenetic trees it is necessary to distinguish between 
15 orthologs and paralogs. The present invention provides a method useful for 

discriminating between orthologs and paralogs and inferring the existence of as yet 
unidentified genes and/or proteins. 

The present invention relates to natural language processing and 
extraction of relational information associated with genes and proteins that are found 
20 in genomics journal articles. Specifically, the natural language processing system of 
the invention is used to parse the articles published in biological journals focusing on 
structure and interactions among genes and proteins followed by computer 
representation of such interactions. 

In accordance with the present invention, specialized databases are 
25 developed that contain information on gene/protein structure and interactions based on 
information derived from preexisting databases and/or research articles including 
information on interactions among genes and proteins, their domain/motif structure 
and their subcellular and tissue expression/distribution patterns. 

The invention relates to a sequence analysis program which utilizes the 
30 specialized database for comparison of a single sequence, processing the output into a 
sequence alignment, computing phylogenetic trees, and analyzing these trees to 
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predict undiscovered genes. This program also includes a set of tools for generating 
motif/domain models from multiple sequence alignments of known genes and for 
using these models for extraction of structurally and/or functionally homologous 
sequences from databases which contain raw sequence data. 
5 The invention further provides for a simulation and hypothesis testing 

program which relies on the specialized databases of gene/protein interactions for 
identifying potentially undiscovered members of multigene families through 
comparisons of regulatory networks for different species and testing hypotheses with 
regard to regulatory cascades. A comparison of homologous regulatory networks 
10 within the same organism and between different species of organisms will allow the 
identification of genes absent in one of the systems under comparison, thus providing 
a set of candidate genes. In this way, genes that contribute to the phenotype of a 
specific disease associated with a particular biological system under analysis may be 
identified, mapped and subjected to mutational analysis and functional studies. 



15 4. BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a block diagram illustrating the three major programs of the 
method according to the present invention: (i) the generation of specialized databases 
based on information on gene/protein structure, function and regulatory interactions 
derived from research papers and databases; (ii) sequence analysis; and 
20 (iii) simulation and hypothesis testing; 

Figure 2 is a block diagram of an information extraction system in 
accordance with a preferred embodiment of the present invention; 

Figure 3 is a diagram illustrating the object representation of molecules 
and relations between them; 
25 Figure 4 shows a set of keywords defining proteins involved in 

apoptosis pathways, these keywords having been utilized for generating a specialized 
sequence database Apoptosis3, this list having been compiled manually for testing the 
concept of specialized databases; 

Figure 5 shows a "species tree," which is a graph depicting the correct 
30 order of speciation events leading to a set of present day species; a "gene tree," which 
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is a graph depicting a history of a few genes from the same species, where each 
species can be represented by multiple paralogous genes (because the set of known 
genes is incomplete for most genomes, and there are often multiple representations of 
the same gene family in the same genome, the gene tree can be drastically different 
5 from the corresponding species tree); and a "reconciled tree", which is the gene tree 
that would be obtained if gene deletions were completely forbidden and all genes were 
known for all species under analysis; 

Figure 6 shows the original tree of ALDH sequences, indicating 
sequence clusters where bacterial, plant, fungal and nematode orthologous genes are 
1 0 present, but a human ortholog was not yet known; 

Figure 7 shows the same phylogenetic tree as in Figure 6 with an 
additional human protein, referred to as antiquitin which was discovered by the 
method of the invention; 

Figure 8 is a schematic diagram illustrating functional network-based 
1 5 gene discovery in accordance with the present invention; 

Figure 9A presents diagrams depicting the regulatory relationships 
among hypothetical proteins (denoted with Arabic numerals) of hypothetical species 
A and B. Proteins in different species denoted with the same numeral are considered 
orthologous. The diagrams show that regulatory relationships between a pair of 
20 proteins can be of three different kinds; 

Figure 9B, 9C, and 9D are diagrams representing Boolean operations 
OR, AND, and XOR, on arcs of the two oriented graphs of Figure 9A, the same 
operations being applicable to the set of vertices of the two oriented graphs; 

Figure 10 is a diagram representing a hypothetical example of defining 
25 homologous protein networks in two different species using protein motifs, the 

diagram showing only two hypothetical proteins (1 and 2) for species A and three 
hypothetical proteins (1,3, and 4) for species B. Protein 1 in both species has motifs 
a and P, protein 2 has motifs 6, e, and £, and proteins 3 and 4 have motifs 8 and (, 
and e, respectively. The motif analysis can indicate that proteins 3 and 4 in species B 
30 may collectively perform the same function as protein 2 in species A; 
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Figure 1 1 A and 1 IB are diagrams respectively representing 
hypothetical examples of evaluating the impact of a "knockout" of hypothetical gene 
A on the expression of a hypothetical gene B. The effect of knock-out of gene A 
calculated by multiplication along the shortest pathway connecting genes A and B is 
5 inhibition of gene B, the resulting effect being zero if the orientation of only one arc in 
the same pathway is reversed; 

Figure 12 is a flow chart representing the scheme of gene discovery 
analysis involving motif/domain analysis in accordance with the present invention; 
and 

10 Figure 13 Identification of genes in C elegans containing either POZ 

or kelch domains. The protein excession numbers are indicated adjacent to the 
different protein domains, the protein corresponding to accession number gi/1 132541 
contains a POZ domain, death domain, kinase domain and heat repeat. 

Figure 14A. Two human sequences with the closest homology to the C. 
1 5 elegans sequence gi/1 1 3254 1 . 

Figure 14B. Computed gene tree indicating that the identified human 
gene represents an ortholog of the C. elegans gene gi/1 132541 . 

Figure 14C. Nucleotide sequence of the death domain gene. 

Figure 14D. Deduced amino acid sequence of the death domain 

20 protein. 

Figure 15. Identification of candidate gene implicated in the etiology of 
Chronic Lymphocytic Leukemia (CLL). Sequence homology between a CLL region 
open reading frame and mouse Rptl (sp/P15533/RPTl) is presented. 

Figure 16A-B. Model of regulatory functions of Rptl. Figure 16A 
25 indicates that in mouse T lymphocytes Rptl serves as a repressor of the gene for 

interleukin 2 receptor (IL-2R). Figure 16B demonstrates that when Rptl is knocked 
out, the regulatory effect is manifested as a block of the apoptotic pathway for T- 
lymphocytes resulting in accumulation of T-lymphocytes in blood. 

Figure 17A. Two EST sequences identified by searching a protein 
30 dbEST using the mouse Mad3 protein as a query. 
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Figure 17B. Nucleotide sequence of the human Mad3 gene. 

Figure 17C. Complete sequence of the human Mad3 protein. A search 
was conducted to identify overlapping sequences. The complete sequence of the gene 
was assembled and the amino acid sequence deduced. The translated human Mad3 
5 sequence consists of 206 amino acid residues 81% of which are identical to the mouse 
Mad3 protein. 

Figure 17D. Multiple alignment of the human Mad3 amino acid 
sequence with known Mad proteins. 

Figure 18A. Phylogenetic tree indicating relationship between three 
1 0 known mouse Mad genes and their two human homologs. 

Figure 18B. Phylogenetic tree including new human Mad3 sequence. 
The phylogenetic tree indicates that the new human gene belongs to the family of Mad 
proteins and is an ortholog of mouse Mad3. 

15 5. DETAILED DESCRIPTION OF THE INVENTION 

The present invention provides methods for identification of novel 
genes comprising: (i) generating specialized databases containing information on 
gene/protein structure, function and regulatory interactions and, (ii) sequence analysis 
which includes homology searches and motif analysis thereby identifying a putative 

20 novel gene of interest. The invention may further comprise performing simulation and 
hypothesis testing to identify or confirm that the putative gene is a novel gene of 
interest. 

The specialized databases are constructed utilizing information 
concerning gene/protein structure or function derived from unpublished data, research 

25 articles and/or existing databases. The specialized databases can be used to identify 
novel genes by: (i) searching for motif/domain combinations characteristic for a 
putative gene of interest; (ii) phylogenetic tree analysis of homologous genes for 
predicting the existence of yet undiscovered genes; (iii) comparing members of 
interactive gene/protein networks from different species for predicting the existence of 

30 yet undiscovered genes; and (iv) testing a hypothesis with regard to known 
interactions of homologs from other species in regulatory pathways. 
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5.1- THE NATURAL LANGUAGE PROCESSING 
The present invention relates to a natural language processing system 
that is designed to parse the electronic versions of articles published in journals that 
report on structural interactions among genes and proteins. The system provides a 
5 method for extracting information on interactions among genes and proteins, their 
domain/motif structure, and/or their sub-cellular and tissue expression/distribution 
patterns, followed by computer representation of such information. 

The general natural language-processing system of the invention is 
schematically depicted in Figure 2. The collection phase automatically collects 

1 0 articles from appropriate literature, and selects articles that contain relevant 

information using Keyword search techniques. In the next phase, the preprocessor 
standardizes the selected articles so that they consist of tagged ASCII text where the 
tags delineate critical components of the article. The next phase, termed the extraction 
phase, retrieves and classifies biological entities, i.e. 9 as names of proteins, genes and 

15 small molecules. In addition, the relationship extraction phase recovers structural 

relationships between the entities. This phase is followed by a phase which performs 
an analysis of the sequence of events. 

The final phase of the system processes the output extracted from an 
article to remove redundancies, inconsistencies and to incorporate implicit 

20 information before adding the extracted knowledge consisting of biological entities, 
their attributes, conditional constraints, and relationships between them, for 
subsequent use in analysis and hypothesis testing. The information extraction system 
as depicted in Figure 2, referred to herein as "GENIE," is designed for use as a general 
processor within the domain of genomics literature although the system may also be 

25 used in other specialized domains. GENIE is an adaptation of MedLEE developed for 
the medical domain. GENIE uses the same source code as MedLEE but the Lexicons 
and grammar were adapted for genomics literature. 

The information extraction system of the present invention is described 
below, by way of example, with reference to the genomics domain uses of GENIE. It 

30 is written in Quintus Prolog and uses the Unix or Windows operating systems, as 
described in detail below. 



WO 00/63687 



PCT/USOO/10302 



13 

A natural-language phrase included in text document is understood as a 
delimited string comprising natural-language terms or words. The string is computer 
readable as obtained, e.g., from a pre-existing database, a keyboard input, optical 
scanning of typed or handwritten text, or processed voice input. The delimiter may be 
5 a period, a semicolon, an end-of-message signal, a new-paragraph signal, or any other 
suitable symbol recognizable for this purpose. Within the phrase, the terms may be 
separated by another type of delimiter such as a blank or another suitable symbol. 

As a result of phrase parsing, terms in a natural-language phrase are 
classified, (e.g., as referring to a gene, a protein, or their interactions) and the 
10 relationships between the interactions are established and represented in a standard 
form. For example, in the sentence "Rap inhibited fyn", the structured form would 
be: 

[action,inactivate,[protein,rap],[protein,fyn]]. 
In such an example, the interaction is "inactivate", the agent is "Rap" and the target 
15 is "fyn." More complex sentences consisting of nested relationships, such as "The 
activation of BAD was suppressed by the phosphorylation of JNK" can also be parsed 
and represented appropriately. The structured output form for this sentence would be: 
[action,inactivate,[action,phosphorylate,x,[prote^ 

] 

20 In the first example, the primary interaction is "inactivate"; in the second example, an 
interaction "phosphorylate" is the agent where the protein "jnk" is its target (the agent 
of "phosyphorylate" in not specified and thus is represented as "x"). In this example, 
the target of "inactivate" is also an interaction "activate" where the target is the 
protein "bad" and the agent is unknown. 

25 While parsing is based on both syntactic and semantic grammatical 

patterns, the substances in a domain are normally only semantic categories such as 
"protein", "gene", and "small molecule." There are no corresponding syntactic 
categories needed for these substances because they are normally all nouns. However, 
each action can be categorized both semantically and syntactically. An action, which 

30 is a semantic category, can generally occur syntactically as a verb "inactivate" or as a 
noun "inactivation." Therefore there are two sets of lexical entries for the actions: 
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syntactic and semantic. The syntactic lexicon for actions specifies the main syntactic 
category such as "v" for verb, "ving" for progressive form of verb, and "activation" 
for noun. The semantic entries for actions not only categorize the actions, but also 
specify features for each action. For example, one feature provides the number of 
5 arguments that are expected for the action, i.e. , some actions are associated with two 
arguments because they have an agent and a target as "inactivate", and others just 
have an agent "mutate." The lexicon of substances and structures appears as 
Appendix A; the syntactic lexicon for actions appears as Appendix B; and the 
semantic lexicon of actions appears as Appendix C. 

10 A second feature specifies whether or not the arguments should be 

reversed when obtaining the target form. For example the arguments of "attributable 
to" should be reversed, i.e., in "the phosphorylation of jnk is attributable to the 
activation of bad", the underlying action is "cause" (from "attributable to"), the agent 
is the "activation of bad" and the target is "the phoshorylation of jnk"), whereas the 

15 arguments of "activates" is not( i.e. in "jnk activates bad" , the agent is "jnk" and the 
target is "bad"). 

Figure 2 shows a preprocessor module of GENIE by which natural- 
language input text is received. The preprocessor thus performs lexical lookup to 
identify and categorize multi-word and single word phases within each sentence. The 

20 output of this component consists of a list of word elements where each element is 
associated with a word or multi-word phrase in the report. For example, assuming 
that the sentence "bad functions as a negative regulator of the activation of jnk" is at 
the beginning of the report, it would be represented as a list of elements where each 
element is a word or phrase. For example, element 1 is associated with "bad", 

25 element 2 with the multi-word phrase "functions as a negative regulator of, element 8 
with "the", and element 9 with "activation". The remainder of the list of word 
positions would be associated with the remaining words in the report. Some of the 
phrases may not need lexical lookup because they already have been tagged by a 
previous component. Such a tagging system is described below in Section 5.2. 

30 The second component of the GENIE system is the parser. It utilizes 

the grammar and categories assigned to the phrases of a sentence to recognize well- 
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formed syntactic and semantic patterns in the sentence and to generate structured 
output forms. The parser proceeds by starting at the beginning of the sentence element 
list and following the grammar rules. When a semantic or syntactic category is 
reached in the grammar, the lexical item corresponding to the next available 
5 unmatched element is obtained and its corresponding lexical definition is checked to 
see whether or not it matches the grammar category. If it does match, the word or 
phrase is removed from the unmatched sentence list, and the parsing proceeds. If a 
match is not obtained, an alternative grammar rule is tried. If no analysis can be 
obtained, an error recovery procedure is followed so that a partial analysis is 

1 0 attempted. The actual grammar used for GENIE appears as Appendix D. 

The parser module of GENIE uses the lexicon, and a grammar module 
to generate target forms. Thus, in addition to parsing of complete phrases, subphrase 
parsing can be used to an advantage where highest accuracy is not required. In case a 
phrase cannot be parsed in its entirety, one or several attempts can be made to parse a 

15 portion of the phrase for obtaining useful information in spite of a possible loss of 
information. 

Conveniently, each module is software-implemented and stored in 
random-access memory of a suitable computer, e.g., a work-station computer. The 
software can be in the form of executable object code, obtained, e.g., by compiling 

20 from source code. Source code interpretation is not precluded. Source code can be in 
the form of sequence-controlled instructions as in Fortran, Pascal or "C", for example. 
Alternatively, a rule-based system can be used such a Prolog, where suitable 
sequencing is chosen by the system at run-time. 

An illustrative portion of the GENIE system is shown in the Appendix 

25 D in the form of a Prolog source listing with comments. The following is further to 
the comments. 

Process jsents with getjnputsents, processjfects and outputresults 
reads in an input stream, processes sections of the input stream according to parameter 
settings, and produces output according to the settings, respectively. Among 
30 parameters supplied to Process jsents are the following: Mode (specifying the parsing 
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mode) and Protocol (html or plain). Process jsents is called by another predicate, after 
user-specified parameters have been processed. 

The parsing modes are selected by GENIE so as to parse a sentence or 
phrase structure using a grammar that includes one or more patterns of semantic and 
5 syntactic categories that are well-formed. For example, for the phrase "bad inactivates 
jnk" a legitimate pattern can be substance 1 action substance2, wherein substance 1 = 
protein bad, action = "inactivates" and substance2 = "jnk." However, if parsing fails, 
various error recovery modes are utilized in order to achieve robustness. The error 
recovery techniques use methods such as segmenting the sentence, processing large 

10 chunks of the sentence, and processing local phrases. Each recovery technique is 

likely to increase sensitivity but decrease specificity and precision. Sensitivity is the 
performance measure equal to the true positive rate of the natural language processing, 
le, , the ratio of information extracted by the natural language processing system that 
should have been extracted. Specificity is the performance measure equal to the true 

1 5 negative information rate of the system, /. e. , the ratio of information not extracted by 
the NLP system that should not have been extracted. Precision is the reliability of the 
system, i.e., the ratio of information extracted correctly compared to all the 
information that was extracted. In processing a report, the most specific mode is 
attempted first, and successive less specific modes are used only if needed. 

20 In accordance with the preferred embodiments of the present invention, 

the parser of Figure 2 includes five parsing modes, Modes 1 through 5, for parsing 
sentences or phrases. Nominally, the parser is configured to first select Mode 1 . If 
Mode 1 is not possible, the program continues with Mode 2 and so forth until parsing 
is complete. With Mode 1 , the initial segment is the entire sentence and all words in 

25 the segment must be defined. This mode requires a well-formed pattern for the 
complete segment. 

Mode 2 requires that the sentence or phrase be segmented at certain 
types of words or phrases, e.g., ** is attributable to." Here, an attempt is made to 
recognize each segment independently, Le. 9 a first segment ending with the word "is " 
30 and a second segment beginning with the word after "to." The segmenting process is 
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repeated until an analysis of each segment is obtained or until segmenting is no longer 
possible. 

Mode 3 requires a well-formed pattern for the "largest'* prefix of the 
segment, i.e., usually at the beginning of the segment. This occurs when a sentence 
5 contains a pattern at the end which is not in the grammar but a beginning portion that 
is included. For example, in "bad inactivates jnk at this time", the beginning of the 
sentence "bad inactivates jnk" will be parsed and the remainder will be skipped. 

Mode 4 requires that undefined words be skipped and an analysis be 
attempted in accordance with Mode 1 . Mode 4 is useful where there are 
10 typographical errors and unknown words. For example, in the phrase "abc bad 

inactivates jnk", the word abc is unknown to the system and will be ignored but the 
remainder of the phrase will be parsed. 

Mode 5 first requires that the first word or phrase in the segment 
associated with an action be found. Next, an attempt is made to recognize the phrase 
1 5 starting with the leftmost recognizable argument. For example, in "during bad 

inactivates jnk on the fifth day," the phrase "bad inactivates jnk" will be parsed and 
the remaining words will not be. If no analysis is found, recognition is retried at the 
next possible argument to the right. This process continues until an analysis is found. 

Processjsects with get section and parse ^sentences gets each section 
20 and generates intermediate output for the sentences in each section. 

Write produces the output as a list consisting of relations and 

interactions 

Setargs sets arguments or parameter values based on user input or by 

default. 

25 The structured output generated by the GENIE program uses a frame- 

based representation. Each frame specifies the informational type, the value, and 
arguments or modifier slots which are also frames. Consider the text data input "bad 
inactivates the phosphorylation of jnk." A corresponding output, as shown below, is a 
frame denoting an action, which has the value inactivate; in addition, there are two 

30 arguments. The first argument is a protein bad and the second argument is an action 
with the value phosphorylate, which has two arguments. The first argument is x 
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signifying that the agent has not been specified; the second argument is a protein with 
the value jnk. The second argument is the target: 

[action,inacti ve, [protein,bad] , [action 5 phosphory late, x, [protein j nk 
In summary, a computer system has been disclosed that generates 
structured information concerning protein and gene interactions and relationships. 



5.2. USE OF BLAST FOR FINDING GENE AND 
PROTEIN NAMES IN JOURNAL ARTICLES 
In a specific embodiment of the invention, an exhaustive list of gene 
and protein names, extracted from GeneBank, is translated into a different alphabet 

1 0 system by substituting each character in the name with a predetermined unique 
nucleotide combination. The encoded names are then imported into the BLAST 
database using the FASTA format. The scientific journals are translated, using the 
same nucleotide combinations, into a continuous string of nucleotides. A query is then 
used to match the translated journals against the nucleotide representation of gene and 

1 5 protein names in the BLAST database. Significant alignments associated with gene 
and protein names are listed in the BLAST output file, which is subsequently 
processed using Perl-scripts. The final result consists of the original journal article 
with XML tags surrounding the gene and protein names. 

To adapt the problem to BLAST'S statistical foundation, different 

20 measures were undertaken to limit the output to the most relevant gene and protein 
names. In addition, in order to fine-tune the matching process, different BLAST 
parameters were adjusted, such as the word size (which sets the size of the high 
scoring words, thus influencing the sensitivity of finding HSPs) and mismatch penalty 
(exact vs approximate matching). 

25 In a specific embodiment of the invention, gene and protein names are 

extracted from GeneBank's gene symbol index file. The following is an excerpt of the 
file after discarding entries that are either composed of only numbers or of less than 
two alphabetic letters: 

gfap gamma 

30 hoxalO 
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hox al 
wac 3'-end 
pit-l/ghf-1 variant 
[...] 



5 This list of gene and protein names is translated into a different 

alphabet system by substituting each character in the name with a predetermined 
unique nucleotide combination. The conversion chart is listed in Appendix E. The 
encoded names are then imported into the BLAST database using the FASTA format. 
For example, the first entry in the list above is "gfap gamma." After translation using 
10 the conversion chart, the same name appears as follows: 

AGCAACTAAACACCCATCCAAGCAAACACACACACAAAC 
Thus, the complete FASTA entry looks like this: 
>gi|l species,gp,gfap gamma 

AAGCAACTAAACACCCATCCAAGCAAACACACACACAAAC 

1 5 In FASTA, the definition line (marked with f > ? ) contains information 

about the database entry. This line can contain any kind of information. The 
information important for this particular example is the third entry in the definition 
line, f gp', that specifies that the name can represent a gene or a protein. If the name is 
unambigous, then the definition line states that the name is only associated with a 

20 gene ('g') or protein ('p'). The fourth entry in the definition line is the name of the 
protein or gene, "gfap gamma" in this case. 

The second line in the FASTA format normally contains the actual 
sequence of the protein/gene. In the example presented, the second line contains the 
translated protein or gene name. 

25 All gene and protein names are translated into the nucleotide 

representation and converted into the FASTA format. Then, the database containing 
these FASTA entries are specially compiled for use in BLAST queries using a 
program that is included in the BLAST package called "formatdb". 

Thus, the scientific journals are translated, using the same nucleotide 

30 combinations, into a continuous string of nucleotides. For example, the sentence "In 
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the absence of costimulation, T cells activated through their antigen ..." is translated 
into 

"AAGTACAGATCCACGGAAGGAACGATCCAAACAAAGACGCAACGACAG 
AAATAACGATCCACATAACTATCCAAATACATACGCACGGAAGTACACAC 
5 GTAATTAAACACGGAAGTACATACAGATCCATCCACGGATCCAAATAACG 
AATTAATTACGCATCCAAACAAATACGGAAGTACTCAAACACGGAACGAA 
CCATCCACGGAAGGACCTACATACGTAAGCAAGGATCCACGGAAGGAAC 
GAAGTACCTATCCAAACACAGACGGAAGTAAGCAACGACAGATCC H 

A query is then used to match the translated journals against the 
1 0 nucleotide representation of gene and protein names in the BLAST database. The 
query is executed using the blastall program that is included in the BLAST package. 
The query line looks like: 

blastall -p blastn -d FASTA.dat -i query.txt 
The flag 'p' denotes the sub-program (blastn is a sub-program of 
1 5 blastall that performs nucleotide matches), 'd' denotes the file that contains the 
FASTA entries and ? denotes the translated query text. 

Significant alignments associated with gene and protein names are 
listed in the BLAST output file. This is an excerpt from a BLAST output file: 

gi|63624 species,gp,ner 
20 Length = 12 

Score = 24.4 bits (12), Expect = 3e-05 
Identities == 12/12(100%) 
Strand = Plus / Plus 
Query: 729 acagaacgacct 740 
25 Sbjct: 1 acagaacgacct 12 

The first line denotes the database entry. The second line denotes the 
database sequence length, followed by the alignment score and the E- value. The next 
line indicates paired matches, mismatches and gapped alignment (the latter two are 
not shown in this example). The lines 'Query* and 'Sbjcf show the actual alignment 
30 between the query and database sequence. This output file is subsequently processed 
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using a Perl-script (see Appendix F). The script shown in Appendix G scans the 
output file, which is sometimes several megabytes long, for any segments that start at 
position 1 of the database sequence (thus disregarding any segments that are only part 
of the sequence). In addition, the script allows for 10% mismatches between the 
5 aligned sequences for long sequences (as shown in the script of Appendix E), or 0% 
mismatches for short sequences. After scanning the output file, an intermediary file 
that lists the candidate sequences is created: 

tran|365|381|gp|18493 

tran|l|17|gp| 18493 
10 peci|549|565|gp|58106 

il-2|621|637|gp|82396 

il-2|325|341|gp|82396 

gati|193|209|gp|92088 

prod|641|657|gp|52292 
15 rapl|105|121|gp|49898 

spec|545|561|gp|33183 

crip|385|401|gp|l 18905 

crip|21|37|gp|l 18905 

as|161|177|gp|133961 
20 her|65|77|gp|88411 



The intermediary file lists the name of the sequence, followed by the 
starting and end point in the query sequence (corresponds to where the two sequences 
matched), the semantic class of the name (protein, gene or protein/gene). The last 
number is not considered. 

25 The intermediary file is then scanned by another Perl program 

(Appendix G). This program compares the starting end points with the actual text, 
making sure that the matched name is an 'autonomous' entity in the query text. For 
example, while "per" in " per gene" should be recognized as a gene name, "per" in 
"personal" should not be recognized as a gene name. The program recognizes other 

30 characters than the space character delimiting an 'autonomous' gene or protein name. 
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In addition, the script looks for plurals of words. For example, " interleukins" should 
be recognized as a protein name, although only the singular form, "interleukin", is in 
the database. 

The final result consists of the original journal article with XML tags 
5 surrounding the gene and protein names. This is done using the same script as in 
Appendix G: 

blocked <phr sem^'gp'^T cell antigen receptor</phr> (TCR)- and 
<phr sem="gp ,t >CD28</phr>-mediated <phr sem="gp">IL-2</phr> gene 
transcription. Therefore, <phr sem="gp M >Rapl</phr> functions as a negative 
10 regulator of... 

To adapt the problem to BLAST'S statistical foundation, different 
measures were undertaken to limit the output to the most relevant gene and protein 
names. 

BLAST is sensitive to the search space the program works in. Thus, 

15 given a long query sequence and a large sequence database, matches have a lower 
statistical significance because the chances are higher that the matches could have 
occurred by chance alone. In addition, matches with few letters have a lower statistical 
significance than matches with many letters. In order to find all true matches with any 
significance level, some measures were undertaken to address this problem. For 

20 example, (i) the query sequence was divided into 1 0 equal length parts, /. e. , the 

journal article was divided into 10 parts and 10 different queries are run on each part 
separately; (ii) the sequence database (with the gene and protein names) is separated 
into 5 databases, each containing protein/gene names of different length; (iii) gene 
and protein names with less than 3 letters in the database were 'expanded 1 , Le. 9 spaces 

25 were added at the beginning and the end of the name. Doing so, the statistical 

significance of a match containing a short name was higher. A space does not only 
include an empty character. For example, a gene name '^4" could occur in a journal 
article as "kinin 4 (k4) M . It was therefore important to define several characters as 
substitutes for a space character. The alphabet in Appendix E defines the nucleotide 

30 combination ATCC as such a substitute. 
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Working with nucleotides implies that errors involving reading frames 
must be addressed. For example, working with a code of four letters, the nucleotide 
combination ATCTGTCACG could mean ATCT/GTCA or TCTG/TCAC or 
CTGT/CACG . Since the text is translated into a nucleotide combination, only one of 
5 these possibilities is correct. But BLAST can not distinguish between these solutions, 
le.y BLAST would potentially match a database sequence to a wrong reading frame in 
the query sequence, producing many nonsense results that could compromise the 
significance of true results. 

The solution to this problem is a comma-free code. A comma free code 

10 knows only one correct reading frame. BLAST therefore does not produce any 

nonsense results. A comma-free code consists of only one permutation of a nucleotide 
combination. For example, given the nucleotide combination ATCC and its 
permutations CATC, CCAT and TCCA, only ONE of these permutations would be 
included in a comma-free code. The code in Appendix E does represent a comma free 

15 code. Comma-free codes were discussed in the early days of DNA research (Crick et 
al., Proc. Natl. Acad. Sci. 43:416-421). 

In order to fine-tune the matching process, different BLAST 
parameters must be adjusted, for example: word size (which sets the size of the high 
scoring words, thus influencing the sensitivity of finding HS?s);mismatch penalty 

20 (exact vs approximate matching); numbers of alignments to show (true matches of low 
significance can sometimes be at the very end of the BLAST output, therefore many 
alignments have to be shown); and expectation value (which sets the significance 
value for matches in the output file). 



5.3. GENERATION OF SPECIALIZED DATABASES 
25 In accordance with the present invention, specialized databases may be 

developed that contain information derived from unpublished data, publications such 
as research articles, theses, posters, abstracts, etc. and/or databases concerning 
interactions among genes and proteins, their domain/motif structure, and their 
biological functions. 
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For example, but not by way of limitation, a specialized database may 
be prepared as follows. Protein and gene sequences may be provided, for example, by 
the Java program PsiRetrieve which allows for quick retrieval of protein or nucleotide 
sequences from binary BLAST databases by sequence accession number, keyword or 
5 groups of keywords, or species name. In addition, using the program PsiRetriever, 
sequences encoding the proteins of interest may be retrieved from the non-redundant 
(NCBI) database of protein sequences and stored as a FASTA file. The FASTA file is 
then converted into a binary blast database using the program FORMATDB from the 
BLAST suit of programs. 

10 Known motifs/domains for proteins may also be collected using the 

flat file versions of major protein databases, such as SwissProt 
(http://expasy.hcage.ch/sprot) and the non-redundant database of NCBI 
(http://www3.ncbi.nlm.nih.gov). The databases can be downloaded and searched for 
the keywords "motif 1 and "domain" in the feature tables of proteins. In addition, 

1 5 existing databases of motifs and domains, such as BLOCKS 

(http://dupsas.Wei2maim.ac.il^cd^cdparent^/databanksblocks/hfml) and 
pfam(http://www.sanger.ac.uk//software/pfam; http://pfm.wustl.edu), can be 
downloaded (Henikoff et al., 1991, NAR 19:6565-6572). Still further, it is understood 
that any publically available database containing gene/protein sequences may be 

20 utilized to generate the specialized databases for use in the practice of the present 
invention. 

Homologous sequences may be aligned using, for example, the 
CLUSTALW program (Higgins, et al. 1996 Methods in Enzymology 266: 383-402). 
A protein's sequence corresponding to each domain/motif can be identified, saved 

25 and used for building a Hidden Markov Model (HMM) of the domain/motif using a 
HMMER and HMMER2 packages (see, Durbin, R. et al. 1998 in Biological Sequence 
Analysis: Probablistic Models of Proteins and Nucleic Acids). HMMER and 
HMMER2 packages are useful for (i) building HMMs from sets of aligned protein or 
nucleotide sequences, and (ii) comparing the HMMs with sequence databases aimed 

30 at identifying significant similarities of HMMs with database sequences. Both 
nucleotide and protein databases can be used for this purpose. Alternatives to the 
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Hidden Markov Model method for building domain/motif models include neural 
network motif analysis (Wu, C.H. et al., 1996, Comput Appl Biosci 12, 109-18; Hirst, 
J.D., 1991, Protein Eng 4:615-23) and positional weight matrix analysis (Claverie, 
J.M., 1994, Comput Chem 18:287-94; Venezia, D., 1993, Comput Appl Biosci 9:65- 
5 9; Bucher, P. 1996, Comput Chem, 20:3-23; Tatusov, R.L., 1994, Proc Natl Acad Sci 
USA 91:12091-5). 

Once a comprehensive collection of motifs/domains is created, each 
particular protein may be compared against a complete database of HMMs to identify 
known motifs and domains, 
10 The Hidden Markov Model (HMM) is built using the following steps: 

Al . Start with a motif/domain name and a single amino acid sequence 

representing a domain or motif. 
A2. Do PSI-BLAST (BLASTPGP) search with the motif/domain sequence 
against a protein non-redundant database. 
15 A3. Retrieve the sequences identified in the database search from the 

protein sequence database. Exclude low-complexity sequences, short 
or incomplete sequences and sequences with similarity score above a 
selected threshold of PPD value <0.001 
A4. Align the set of sequences with CLUSTALW (or other multiple 
20 sequence alignment program). 

A5. Use the set of aligned sequences for building HMM with the programs 
provided with HMMER and HMMER2 packages (see Hughey and 
Krogh 1996, J. Mol. Biol. 235:1501-1531). 
A6. Do a new database search comparing new HMM with the non- 
25 redundant protein database. 

A7. Continue steps A3-A6 until the convergence of the Markov model /. e., 
until no new sequences are identified, or the maximum allowed 
number of iterations as defined by the user is reached. (Hugh R. and 
Rrogh A., 1996, Comput. Appl. Biosci. 12: 95-107). 
30 In addition, in yet another embodiment of the invention, a specialized 

database may be designed to contain a semantic model of proteins and of the possible 
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interactions between them. Such databases are particularly useful for computation 
and analysis of regulatory networks between proteins. The semantic model is designed 
for representing substances, such as proteins and actions between them, and is based 
on widely accepted principles of object-oriented programming languages such as Java. 
5 Figure 3 is a diagram illustrating the object representation of molecules and relations 
between them. As indicated in Figure 3 there are six major classes, corresponding to 
the top-level classification of objects and actions: (/) a substance; (ii) a state of a 
substance; (iii) a similarity between substances; (iv) an action between substances; (v) 
a result of the action; and (vi) a mechanism that enables an action. 

10 Figure 3 presents the class design graphically, listing the variables that 

represent the properties of each class or class object in the implementation. Classes 
can be made nested via the mechanism of "inheritance", i.e., classes are defined 
starting with the most general ones and moving towards more specific classes. 
Definition of more specific classes is simplified because the properties of the general 

15 classes are "inherited" by the specific classes and need not be redefined each time 

(see, Flanagan 1997, Java in a Nutshell, Second Edition. O'Reilley & Associates, Inc. 
Sebastopol, CA). 

As shown in Figure 3, the two key object types in this scheme are 
substances (nodes of the graph representing regulatory networks) and actions 

20 (oriented edges connecting pairs of nodes), while result and mechanism objects are 
auxiliary to object action. Each substance object is characterized with a state. In this 
scheme, action is the most complicated object; each action object is characterized by a 
specific pair of substances participating in the action, one of which can be active and 
is referred to as Subject Substance and the second of which can serve as a substrate for 

25 the former and is referred to as Object Substance. Furthermore, for each action the 
initial and final states corresponding to interacting substances are defined. The 
property Time Required of each Action Object allows the setting of different 
durations for different actions (time is measured in relative units; see Rene Thomas 
and Richard D'Ari, 1990, "Biological Feedback," CRC Press Boca Raton, Ann Arbor, 

30 Boston). 
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Once developed, the specialized databases can be used to identify 
novel genes based on computation and analysis of phylogenetic trees for multigene 
families and analysis of homologous regulatory networks. 

In a specific embodiment of the invention, a specialized database was 
5 generated using a set of keywords defining proteins involved in apoptosis (see, Figure 
4). The specialized sequence database was referred to as Apoptosis 3. As a first step 
in generating the specialized database, a comprehensive set of articles describing the 
system of apoptosis or programed cell death was compiled. The articles were 
analyzed and information on regulatory pathways characterizing apoptosis from a 

10 variety of different organisms was extracted. Such pathways included those involved 
in MHC-T cell receptor interactions, inflammatory cytokine signal transduction, 
induction by light, y-radiation, hyperosmolarity or heat shock, pathways involving 
immunoregulatory receptors or receptors having cytoplasmic domains, integrin- 
related pathways and perforin/granzymep related pathways. The collected 

1 5 information was stored using Powerpoint (Microsoft) as a collection of graph/plots 
depicting the regulatory pathway. In addition, a list of proteins relevant to regulation 
of apoptosis was compiled. 

Using the program Psi Retriever, sequences encoding the proteins 
relevant to regulation of apoptosis were retrieved from the non-redundant (NCBI) 

20 database of protein sequences and stored as a FASTA file. The FASTA file was then 
converted to a binary blast database using the program FORMATDB from the 
BLAST suit of programs. The BLAST suit of programs provides a set of programs 
for very fast comparisons of a single sequence to a large database. Both the database 
and the search or query sequence can be any combination of nucleotide and/or amino 

25 acid sequences. 

In a working example described herein, the Apoptosis 3 database was. 
used to compare genomic and cDNA sequences derived from the 13q region of human 
chromosome 13. This region of the chromosome is associated with Chronic 
Lymphocytic Leukemia (CLL). Using this method of analysis a human gene with 

30 significant homology to the mouse Rptl gene was identified. When the activity of 
Rptl is knocked out in mice, the regulatory effect is manifested as a block in T- 
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lymphocyte apoptosis. This result indicates that the identified human Rptl homology 
may represent the gene in which genetic defects lead to CLL. 

The amino acid sequence of the human Rptl gene is presented in 
Figure 15. The present invention relates to nucleic acid molecules encoding the 
5 human Rptl protein shown in Figure 15. The invention also relates to nucleic acid 
molecules capable of hybridizing to a nucleic acid molecule encoding the human Rptl 
protein presented in Figure 1 5 under conditions of high stringency. By way of 
example and not limitation, procedures using such conditions of high stringency are as 
follows: Prehybridization of filters containing DNA is carried out for 8 hours to 

1 0 overnight at 65°C in buffer composed of 6x SSC, 50 mM Tris-HCl (pH7.5), ImM 
EDTA, 0.02% PVP, 0.02% Ficoll, 0.02% BSA and 500 mg/ml denatured salmon 
sperm DNA. Filters are hybridized for 48 h at 65°C in prehybridization mixture 
containing lOOmg/ml denatured salmon sperm DNA and 5-20 x 10 6 CpM of 32 P- 
labeled probe. Washing of filters is done at 37°C for 1 h in a solution containing 2x 

15 SSC, 0.01% PVP, 0.01% Ficoll and 0.01% BSA. This is followed by a wash in 0.1 x 
SSC at 50°C for 45 minutes before autoradiography. Other conditions of high 
stringency which may be used are well known in the art. 

5.4. GENE DISCOVERY THROUGH PHYLOGENETIC 

ANALYSIS OF GENE FAMILIES 

20 The present invention provides a method for identifying novel genes 

comprising the following steps: (i) comparing a single sequence with a database; (ii) 
processing the output into a sequence alignment; (iii) computing gene trees; and (iv) 
analyzing the trees to predict the existence of undiscovered genes. 

Figure 5 shows a "species tree," a "gene tree" and a "reconciled tree". 
25 A "species tree", as defined herein, is a graph depicting the correct order of speciation 
events leading to a set of present day species as defined by taxonomy. A "gene tree" 
is a graphical representation of the evolution of a gene from a single ancestral 
sequence in a common progenitor to a set of present-day sequences in different 
species. Where gene duplication has occurred, a branch is bifurcated. The branch 
30 lengths of a gene tree are most frequently measured either in terms of the number of 
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amino acid or nucleotide replacements per site or in terms of millions of years 
(absolute geological time). In the former case, the average replacement rate in the 
majority of the published trees varies among tree branches, and the root-to-tip 
distances are different for different present day sequences. In the latter case, all root- 
5 to-tip distances are equal and the height of each interior node of the tree corresponds 
to the absolute geological time passed since the gene duplication corresponding to the 
interior node took place. 

If a gene is unique, i.e., represented with a single copy per genome 
rather than being a member of a family of similar genes, the correct gene tree 

10 depicting the origin of this gene in a few different species is identical to the species 
tree. In many instances, a single ancestral gene has been duplicated repeatedly during 
evolution to form a multigene family. A gene tree is constructed from a gene as it 
occurs in several species and reflects both speciation events and gene duplications 
within the same genome. Two homologous genes taken from different species that 

1 5 originated from the nearest common ancestor by speciation are referred to as 

orthologs, while any two genes that originated from the common ancestor via a series 
of events involving intragenomic duplications, or conversions, are called paralogs. 
The terms "ortholog" and "paralog" are applied to both nucleic acid and proteins 
herein. 

20 If gene deletions are forbidden and all genes for all species represented 

in the tree are known, the gene tree can be reconfigured to recapitulate the species 
tree, such that each subtree contains only orthologous genes. This tree is referred to as 
a reconciled tree and is shown in Figure 5. Imperfect gene trees which contain 
incorrect or partial species subtrees can be used to build reconciled trees that indicate 

25 events of speciation, gene loss, and gene duplication. 

Orthologs from different species in gene trees are usually clustered 
together, so that if all the existing homologous genes from different species were 
known, the same relationship of species would be recapitulated in each cluster of 
orthologous genes. Since in reality a considerable number of genes are not yet 

30 identified, the real gene trees contain incomplete clusters of orthologs that can be used 
for identification of the missing genes. 
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By applying phylogenetic analysis, i.e., reconstruction of gene trees of 
gene/protein sequences, one can predict the existence of undiscovered genes in 
humans and other species in addition to identifying the function of a gene. Such a 
technique is a significantly more powerful tool for identification of new genes than 
5 mere sequence comparisons. 

Methods of computing gene trees from a set of aligned sequences 
include the : (i) heuristic method based on an optimization principle which is not 
directly motivated by a probability model (Fitch, 1974 J. Mol. Evol. 3:263-268)), (ii) 
the maximum likelihood method (Goldman, 1990, Syst Zool. 30:345-361; Yang et 
10 al., 1995, Syst. Biol. 44:384-399; Felsenstein, J., 1996, Methods EnzymoL 266-418- 
427); and (iii) the distance matrix tree making method (Saito, N. and Nei, M., 1987, 
Mol. Biol. Evol. 4:406-425). Since the data analyses of orthologs and paralogs often 
involve very distantly related sequences, the maximum likelihood method is 

■ 

preferably used for small data sets and the distance-matrix method in other instances. 

1 5 To construct a reconciled tree according to the invention, the first step 

comprises a search for homologs in a publicly or privately available database such as, 
for example, GenBank, Incyte, binary BLAST databases, Swiss Prot andNCBI 
databases. Following the identification of homologous sequences a global alignment 
is performed using, for example, the CLUSTALW program. From the sequence 

20 alignment a gene tree is constructed using, for example, the computer program 
CLUSTLAW which utilizes the neighbor-joining method of Saito and Nei (1997, 
Mol. Biol. Evol. 4:406-425). Construction of a species tree is then retrieved from, for 
example, the following web site: 
http://www.3.NCBI.NLM.NIH.GOV//taxomy.tax.html. 

25 The species tree and gene tree are given as input into the algorithm 

described below, which integrates both trees into a reconciled tree. Agreement 
between the gene tree and the corresponding species tree for any given set of 
sequences indicates the identification of orthologs. In contrast, disagreement between 
the species and gene tree suggest a gene duplication that resulted in the formation of a 

30 paralog. Thus, through generation of a reconciled tree one can identify orthologs 
present in one species but missing in another. These can be deduced by forming 
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subtrees of orthologs in a gene tree, and then comparing the subtree in the gene tree 
with a species tree. A missing gene appears as a branch present in the species tree but 
absent in the gene tree. The algorithm for defining an orthologous gene subtree and 
predicting the undiscovered, or lost in evolution, genes is as follows: 
5 Let T g be the most likely gene tree identified with one of consistent 

tree-making methods from a set of properly aligned homologous genes {1,2, s}, 
such that one or more homologous genes from every species corresponds to pending 
vertices of T g , Each gene is labeled with the species it comes from (l,...,s) adding 
subscripts to distinguish homologous genes from the same species whenever it is 

1 0 necessary. Let T g be the true species tree (tree correctly reflecting speciation events 
which we assume to be known) for species {1,2, s}. Due to the biological 
meaning of T s each species in this tree is represented only once. It is assumed that 
both T s and T g are binary, although it is straightforward to extend the algorithm 
described here to the case of multifurcated trees. 

15 Algorithm 

Al . For each pair of interior nodes from trees T 8 and T s , compute similarity 
A2. Find the maximum o(S JS 

A3. Save S gi as a new subtree of orthologs, save {S gi } - {S sj } as a set of species 
20 that are likely to have gene of this kind (or lost it in evolution). 

A4. Eliminate S gi from T g ; T g : - T g \S gi . 
A5 . Continue A2 - A4 until T g is non-empty . 
The following definitions apply: 

Let S gi be an ith subtree of T g (corresponding to the ith interior node), 

25 correspondingly, let S sj be yth subtree of tree T r 

Let {S gi } stand for an unordered set of species represented in S gi such that each 
species is represented exactly once, and let \{SJ\ and {\S .1} be the number of 
entries in {S gI } and the number of pending vertices in S gi , respectively. Define by 
S sj (S gi ) the unique subtree of S sj that has leaves labeled exclusively with species from 

30 i ( s gi} I > so that each element of | {S gi } | is used i.e., that is, the unique subtree 
obtained by eliminating from S sj all species that are not present in | {S .} | . 
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Then define similarity measure, a, between S gi and S sj in the following way: 

The support of tree clusters by data can be measured using the bootstrap technique 
5 described in Felsenstein (1985, Evolution 39:783-791). 

In an embodiment of the invention, the human antiquitin gene was 
identified using phylogenetic analysis. The aldehyde dehydrogenase gene family in 
humans can be subdivided into at least ten ancient subtrees characterized by different 
functions of corresponding proteins. These genes probably arose from a series of gene 

1 0 duplications of an ancestral gene which took place before the divergence of a common 
ancestor of Eukaryotes and Eubacteria. 

The aldehyde dehydrogenase gene cluster is highlighted in Figure 6 
which shows the original tree of ALDH sequences, the circled area indicating a 
sequence cluster where bacterial (Bacillus subtilis), plant (Brassica napus\ and 

15 nematode (Caenorhabditis elegans) ortholog is present, but a human ortholog is not 
known. A random screening of cDNA libraries showed that a human ortholog, 
referred to as antiquitin, does exist. Figure 7 shows the same gene tree as in Figure 6 
with an additional human protein referred to as antiquitin present in the tree. 

In yet another embodiment of the invention, a human ortholog of the 

20 murine Max-interacting transcriptional repressor Mad3 was identified through 

phylogenetic analysis of a gene family. The gene tree was constructed as follows. 
The protein sequences of known members of the Mad gene family were extracted 
from GenBank database. The extracted sequences were aligned using multiple 
alignment program CLUSTALW running on Sun SPARC station. Redundant and 

25 non-homologous sequences as well as distant homologs from S> cerevisiae, C. 

elegans, D. melanogaster etc. were removed from the alignment. The refined set of 
sequences were realigned with CLUSTALW and a gene tree as presented in Figure 
1 8 A was computed. To identify a human ortholog of the Mad3 protein, a human 
dbEST at NCBI was searched with program TBLASTN using mouse Mad3 protein 

30 sequences as a query. Two highly homologous ESTs were identified and are 
presented in Figure 17A. To obtain a complete coding sequence a search was 
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conducted to obtain overlapping sequences in dbEST. The search for overlapping 
sequences was performed using the program Iterate with EST Zs77e55.rl 
(gb/AA278224) as the search query. The search identified a single overlapping 
sequence. The search for overlapping sequences was performed using program Iterate 
5 with EST zs77e55.rl (gb/AA278224) serving as a query. The search returned a single 
overlapping sequence, namely HUMGS0012279 (dbj/C02407), thus showing that the 
two EST sequences found during the initial TBLASTIN search belong to the same 
gene. The complete sequence of the gene was assembled from the two ESTs using 
commercially available sequence assembly program SeqMan 1 1 (DNASTAR Inc., WI). 

10 The nucleotide sequence of the human Mad3 gene is presented in Figure 17B. The 
deduced amino acid sequence of which is presented in Figure 17C. The complete 
DNA sequence is also shown. 

The present invention relates to nucleic acid molecules encoding the 
human Mad3 protein shown in Figure 17C. The invention also relates to nucleic acid 

15 molecules that hybridize to the nucleic acid molecule of Figure 17B under conditions 
of high stringency and encode a Mad3 protein. By way of example and not limitation, 
procedures using such conditions of high stringency are as follows: Prehybridization 
of filters containing DNA is carried out for 8 hours to overnight at 65°C in buffer 
composed of 6x SSC, 50mM Tris-HCl (pH7.5), ImM EDTA, 0.02% PVP, 0.02% 

20 Ficoll, 0.02% BSA and 500 mg/ml denatured salmon sperm DNA. Filters are 

hybridized for 48 hours at 65°C in prehybridization mixture containing 100 mg/ml 
denatured salmon sperm DNA and 5-20 x 10 6 CpM of 32 P-labeled probe. Washing of 
filters is done at 37°C for 1 hour in a solution containing 2x SSC, 0.01% PVP, 0.01% 
Ficoll and 0.01% BSA. This is followed by a wash in O.lx SSC at 50°C for 45 

25 minutes before autoradiography. Other conditions of high stringency which may be 
used are well known in the art. 

5.5. SIMULATION AND HYPOTHESIS TESTING 
The simulation and hypothesis testing methods of the invention, 
described in the subsections below, utilize specialized databases of gene/protein 
30 structures and interactions for identifying potentially undiscovered members of 
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multigene families through comparisons of regulatory networks for different species, 
searching expressed sequence tag (EST) databases, and simulation of regulatory 
cascades. 

5.5.1. GENE DISCOVERY THROUGH ANALYSIS 

5 OF REGULATORY NETWORKS 

The present invention provides a method for identifying undiscovered 
genes through comparisons of regulatory networks for different species where 
functionally similar regulatory systems are conserved. The amount of information 
available concerning regulatory genes and/or proteins in different organisms and their 

10 functional relationships allows one to reconstruct and compare regulatory networks. 
Since in most cases, the knowledge of all genes involved in almost any particular 
regulatory system is incomplete, a comparison of homologous networks within the 
same organism and between different species permits the identification of genes 
absent in a system under comparison. 

1 5 The identified genes, being part of a regulatory network, are implicated 

as potentially contributing to a phenotype of a disease associated with the system 
under analysis. Using the methods of the present invention these putative disease 
genes can be cloned, mapped and analyzed for mutations directly, thereby omitting 
the expensive and time-consuming steps of positional cloning and sequencing of 

20 genomic regions. Gene discovery by analysis of regulatory networks is outlined in 
Figure 8. The analysis is initiated starting with a biological system (e.g., signaling 
pathway of genes involved in Bcl-2-regulated apoptosis in lymphocytes), a single 
gene (e.g., Bcl-2) or a gene family (e.g., caspases). 

Initially, a specialized database is generated for comparison of 

25 regulatory networks between different species. For example, starting with a single 
candidate gene in a single species, a typical iteration in this process begins with 
identification of all known proteins and genes that are upstream and downstream with 
respect to it in regulatory hierarchies and the reconstruction of a network of 
interacting genes and proteins. Next, for each protein, a set of key domains and motifs 

30 is identified and this information is used to search for related proteins in humans and 
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other species. The identified sequences are compared and for each pair of sequences 
showing similarity above a certain threshold, a similarity object is generated. A 
similarity object is generated if two sequences, nucleotide or amino acid, show 
significant similarity in database searches (p value < 0.001). The object retains the 
5 following information: (i) reference to similar substances i.e., genes or proteins; 

(ii) significance of the similarity, similarity score and percent of identity; and 

(iii) coordinates of the similarity region within two compared sequences. 
"Orthology objects" constitute a subset of "similarity objects" which satisfies one 
additional requirement, that two similar sequences should be identified as 

10 orthologs by the tree-based algorithm described above. In identifying orthologs, if 
gene A is orthologous to gene B, and gene B is orthologous to gene C, gene A is 
necessarily orthologous to gene C. 

In a specific embodiment of the invention, for each species under 
analysis, orthologous proteins or genes are identified. In a further embodiment of the 

1 5 invention, small orthologous molecules participating in a regulatory network for two 
or more species may also be identified. Where proteins, genes, or molecules are 
orthologs, the action of the protein, gene or molecule between species may be 
interchangeable. If more than two species are involved in the analysis, subtrees of 
orthologous substances and subtrees of orthologous actions are identified. 

20 Once orthologous genes, proteins or molecules are identified in two or 

more species, by forming a reconciled tree, for example, a set of orthologous or 
paralogous regulatory networks can be analyzed and visualized using graph theory 
where arcs represent actions and vertices represent substances. Thus, the method of 
the invention may further comprise the following steps: (i) superimposing the 

25 orthologous regulatory networks from two or more species and searching for the 
actions (arcs) and substances (vertices) in the homologous networks that are 
represented in some taxa but absent in others; (ii) superimposing paralogous 
regulatory networks from the same taxa and searching for paralogous genes that are 
missing in some taxa; and (iii) computing a general regulatory network that 

30 summarizes common regulatory sequence relationships known for more than one 
species. 
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In a specific embodiment of the invention a set of regulatory networks 
from different species, relating to the same biological system, apoptosis, for example, 
can be analyzed and visualized utilizing the following methods: (i) for each species 
functional information is collected relating to apoptosis; (ii) using the functional 
5 information, regulatory networks for each species comprised of interacting proteins 
and/or the genes involved in apoptosis are generated; (iii) the sequences of the 
interacting proteins and genes of each of the regulatory network are compared and for 
sequences showing similarity above a predetermined threshold range; and (iv) 
distinguishing between orthologs and paralogs using the methods set forth above* 

10 An analysis similar to that performed using subtrees of sequences may 

be applied to classify protein functions as orthologous or paralogous actions. A 
"generalized" regulatory network maybe represented as a network wherein a substance 
as it occurs in a particular species is substituted with a cluster (i.e., subtree) of 
orthologous substances among species. In the final step of the analysis the clusters 

15 within each species are compared to one another, to identify missing genes. 

Figure 1 1 depicts the regulatory relationships among hypothetical 
proteins (denoted with Arabic numerals) of hypothetical species A and B. As 
indicated in Figure 1 1 A, an overlay of regulatory data for two species overlaps, but 
not completely. As indicated, protein 5 is known only for species B while protein 3 is 

20 known only for species A. The proteins in different species denoted with the same 
numeral are considered orthologous. As indicated, the regulatory relationships 
between a pair of proteins can be of three different kinds. Figure 9B, 9C, and 9D 
represent Boolean operations, OR, AND, and XOR, as arcs of the two regulatory 
relationships depicted in Figure 9A, the same operations being applicable to the set of 

25 vertices of the two regulatory relationships. 

In some instances, orthologous networks in two distantly related taxa may have the 
same domains but arrangement of the domains between the related taxa may be 
different. In such a case, a one-to-one correspondence between orthologous proteins 
in closely related species has to be substituted with a one-to-many relationship among 

30 domains o comprised within the proteins. For this purpose, a similarity object may be 
defined operating on pairs of motifs/domains in two proteins, and substitute pairs of 
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orthologous proteins with pairs of orthologous domains, After this correction, 
homologous networks are compared as described above. 

Figure 10 is a diagram representing a hypothetical example of defining 
homologous protein networks in two different species using protein motifs, the 
5 diagram showing only two hypothetical proteins (lane 2) for species A and three 
hypothetical proteins (lanes 1, 3, and 4) for species B. Protein 1 in both species has 
motifs a and P, protein 2 has motifs 8, e, and C, and proteins 3 and 4 have motifs 5 
and C, and e, respectively. The motif analysis indicates that proteins 3 and 4 in 
species B may collectively perform the same function as protein 2 in species A. 

10 5.5.2 GENE DISCOVERY BASED ON PROTEIN 

MOTIF/DOMAIN SEARCHES 

The present invention provides yet another method for identifying 
genes that are homologous and perform the same or an analogous function in different 
species. The method of the invention comprises the following steps: (i) creating a 

1 5 database of sequences which comprise a motif or domain composition of a gene of 

interest using, for example, HMMER software; and (ii) searching additional databases 
for expressed sequence tags (ESTs) containing the domains and motifs characteristic 
for the gene of interest with HMMs of domains and motifs identified in step (i). In yet 
another embodiment of the invention, sequences may be searched which correspond 

20 to nucleotide sequences in an EST database or other cDNA databases using a program 
such as BLAST and retrieving the identified sequences. In an optional step, for each 
EST identified, sequence databases can be searched for overlapping sequences for the 
purpose of assembling longer overlapping stretches of DNA. Once identified, the 
ESTs can be used to isolate full length nucleotide sequences comprising the gene of 

25 interest using methods such as those described in Section 5.4, infra. 

The general flowchart scheme for gene discovery analysis based on 
motif/domain search is shown in Figure 11. In a specific embodiment of the 
invention, the method referred to as the "phylogenetic reflection technique"comprises, 
first, defining the motif or domain composition of a gene of interest involved in a 

30 biological system of interest. Second, protein-coding genes from other species, 
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including for example yeast and/or nematode genes, that bear a significant similarity 
to the gene of interest or a specified domain of the corresponding protein are 
collected. Third, the identified genes are in turn subjected to a "domain analysis" to 
establish protein motifs which might suggest a function of these genes using, for 
5 example, HMMER software. Fourth, the selected genes are in turn used for database 
searches in EST databases (dbEST) and/or a non-redundant (nr) database to identify 
unknown genes that are potentially orthologous to the selected yeast and nematode 
genes. Once identified ESTs having different tumor suppressor domains may be 
linked using multiple PCR primers. Using routine cloning techniques, well known to 
10 those of skill in the art, a full length cDNA representing the gene of interest can be 
obtained. 

Once new genes are identified by domain/motif analysis experimental 
searches may be carried out to isolate complete coding sequences and evaluate their 
tissue- and disease-specific expression patterns. In parallel their position with respect 

1 5 to regulatory networks can be identified as described below. 

In a specific embodiment of the invention, an apoptosis related human 
gene was identified using the method described above. As a first step C. elegans 
genes containing either POZ or iCelch domains were identified. A Hidden Markov 
Model was developed using POZ and Kelch sequences from the Drosophila Kelch 

20 protein and any identified homologs. The resulting Hidden Marker Model was used 
to search through the collection of C. elegans protein sequences. One of the identified 
C elegans genes contained a POZ domain, death domain, kinase domain and heat 
repeat. The presence of both a death domain and a kinase domain suggested that the 
protein functions as a regulatory protein. 

25 A human EST database was searched using the protein sequence of the 

identified C. elegans gene and two sequences were identified (Figure 14A). A gene 
tree was computed to determine whether the identified human sequences were 
orthologs of the C. elegans gene. As depicted in Figure 14B, the human EST 
AA481214 appears to be a true ortholog of the C. elegans gene. Figure 14C presents 

30 the nucleotide sequence of the identified death domain gene. Figure 14D presents the 
amino acid sequence of the death domain protein. 



WO 00/63687 



PCT/US00/10302 



39 

The present invention encompasses the nucleic acid molecule of Figure 
14C, comprising the sequence of EST AA481214 and proteins encoded by said 
nucleic acid molecule. The invention also relates to nucleic acid molecules capable of 
hybridizing to such a nucleic acid molecule under conditions of high stringency. By 
5 way of example and not limitation, procedures using such conditions of high 

stringency are as follows: Prehybridization of filters containing DNA is carried out 
for 8 hours to overnight at 65°C in buffer composed of 6x SSC, 50mM Tris-HCl 
(pH7.5), ImM EDTA, 0.02% PVP, 0.02% Ficoll, 0.02% BSA and 500 mg/ml 
denatured salmon sperm DNA. Filters are hybridized for 48 hours at 65°C in 
10 prehybridization mixture containing 100 mg/ml denatured salmon sperm DNA and 5- 
20 x 10 6 CpM of 32 P-labeled probe. Washing of filters is done at 37°C for 1 hour in a 
solution containing 2x SSC, 0.01% PVP, 0.01% Ficoll and 0.01% BSA. This is 
followed by a wash in O.lx SSC at 50°C for 45 minutes before autoradiography. 
Other conditions of high stringency which may be used are well known in the art. 

15 

5.5.3. SIMULATION OF REGULATORY CASCADES 
In an embodiment of the invention, an interactive graphical program is 
utilized for visualizing the scheme of regulatory relationships, "current" states of the 
substances, and active and inactive actions between pairs of substances. Such a 

20 program can be utilized for identification of genes which are associated with a specific 
disease. Currently, disease associated genes are discovered through positional cloning 
methods which combine methods of genetics and physical mapping with mutational 
analysis. The present invention provides a novel method for discovering disease 
associated genes. For simulating regulatory cascades, it is assumed that the time in a 

25 simulated regulatory system advances in discrete "quanta," or periods of time. The 
"state of substances" of the system for each discrete period of time is computed by: 
creating a set of substance objects, where a set of interactions between each created 
substance object is known, an initial state is specified. The time is initially set to zero. 
All defined actions are observed to confirm that the substances corresponding to the 

30 actions (i) exist, and (ii) are in the right initial states. Action is defined by a pair of 
substances that are in suitable states. The "subject" substance is in the inactive state, 
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while the "object" substance can be in either active, or inactive, state depending on the 
action type. For example, the action "dephosphorylation" requires an active 
phosphatase ("subject" substance) and a phosphorylated substitute protein ("object" 
substance) in phosphorylated form. If both conditions are satisfied, the action is 
5 recorded as in progress. At termination, the substances must change their states as 
specified by the action. On each following "quantum" of time, the simulation 

to 

proceeds in the same way while maintaining the "bookkeeping" of the remaining time 
for each action and the remaining lifespan of each substance. The simulation stops 
when there are no more active actions available. The program allows editing of the 

10 properties of the objects, changing the scale and focus of the visualized simulation, 
and experimenting with the systems output. 

In a specific embodiment of the invention a "knock out" of a gene can 
be simulated to model the regulatory system that normally includes hypothetical 
gene A. One of the typical questions related to the gene knock out is how does the 

1 5 knock out affect a biological pathway of interest. A hypothetical example of 

evaluating the impact of a knock out of hypothetical gene A on the expression of a 
hypothetical gene B is shown in Figure 12. The answer to such a question could be 
"gene B will be inhibited" or "gene B will be induced" or "no effect". 

In the practice of the present invention, a simple algorithm involving 

20 multiplication of gene interaction "signs" along the shortest pathway between the 
genes can be used to determine the outcome. The algorithm involves the following 
steps: (i) identification of the shortest non-oriented pathway connecting genes A and 
B involved in a pathway of interest; (ii) assigning sign "-" to gene A since it is 
knocked out and taking this sign as the initial sign value; (iii) moving along the 

25 shortest pathway between genes A and B, multiplying the current value of the sign 
with the sign of the next arc, where stands for inhibition, "+" stands for induction 
or activation, and "0" stands for the lack of interaction between two proteins in the 
specified direction; (iv) determining if the final result of multiplication is "0", if so 
eliminating the zero arc and trying to find the shortest oriented bypass pathway 

30 between A and B in the remaining network; otherwise stop. The final value of the 

sign at the moment of arriving at vertex B would indicate the most likely effect of the 
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knock out of gene A which can be any one of the following: inhibition of gene B, 
induction/activation of gene B, or none. In addition to the "electronic knock out", an 
"electronic knock in" of a particular gene can be simulated. In such a computer 
simulation, the artificial addition of a gene and its effect on a regulatory system may 
5 be analyzed. 



5.6. IDENTIFICATION AND ISOLATION OF NOVEL GENES 
The present invention relates to identification of novel genes, i.e., 
missing orthologs or paralogs, and the isolation of nucleic acid molecules encoding 

10 novel genes. In a specific embodiment, a nucleic acid molecule encoding a missing 
ortholog or paralog can be isolated using procedures well known to those skilled in 
the art (See, for example, Sambrook et al., 1989, Molecular Cloning, A Laboratory 
Manual, 2d Ed., Cold Spring Harbor Laboratory Press, Cold Spring Harbor, New 
York Glover, D.M. (ed.), 1985, DNA Cloning: A Practical Approach MRL Press, 

1 5 Ltd., Oxford, U.K. Vol. I, II.). 

For example, genomic and/or cDNA libraries may be screened with 
labeled DNA fragments derived from a known ortholog or paralog from a specific 
species and hybridized to the genomic or cDNA libraries generated from a different 
species. For cross species hybridization, low stringency conditions are preferred. For 

20 same species hybridization, moderately stringent conditions are preferred. Any 
eukaryotic cell potentially can serve as the nucleic acid source for the molecular 
cloning of the gene of interest. The DNA may be obtained by standard procedures 
known in the art from cloned DNA (e.g., a DNA "library"), by cDNA cloning, or by 
the cloning of genomic DNA, or fragments thereof, purified from the desired cell 

25 By way of example and not limitation, procedures using conditions of 

low stringency are as follows (see also Shilo and Weinberg, 1981, Proc. Natl. Acad. 
Sci. USA 78:6789-6792; and Sambrook et al. 1989, Molecular Cloning, A Laboratory 
Manual, 2d Ed., Cold Spring Harbor Laboratory Press, Cold Spring harbor, New 
York): Filters containing DNA are pretreated for 6 h at 40 °C in a solution containing 

30 35% formamide, 5X SSC, 50 mM Tris-HCl (pH 7.5), 5 mM EDTA, 0.1% PVP, 0.1% 
Ficoll, 1% BSA, and 500 mg/ml denatured salmon sperm DNA. Hybridizations are 
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carried out in the same solution with the following modifications: 0.02% PVP, 0.02% 
Ficoll, 0.2% BSA, 100 mg/ml salmon sperm DNA, 10% (wt/vol) dextran sulfate, and 
5-20 X 10 6 cpm 32 P-labeled probe is used. Filters are incubated in hybridization 
mixture for 18-20 h at 40°C, and then washed for 1.5 h at 55 °C in a solution 
5 containing 2X SSC, 25 mM Tris-HCl (pH 7.4), 5 mM EDTA, and 0.1% SDS. The 
wash solution is replaced with fresh solution and incubated an additional 1.5 h at 
60°C. Filters are blotted dry and exposed for autoradiography. If necessary, filters 
are washed for a third time at 65-68 °C and reexposed to film. Other conditions of 
low stringency which may be used are well known in the art (e.g., as employed for 

10 cross species hybridizations). 

In another specific embodiment, a nucleic acid which is hybridizable to 
a nucleic acid under conditions of moderate stringency is provided. For example, but 
not by way of limitation, procedures using such conditions of moderate stringency are 
as follows: filters containing DNA are pretreated for 6 h at 55° C in a solution 

15 containing 6X SSC, 5X Denhaifs solution, 0.5% SDS and 100 mg/ml denatured 

salmon sperm DNA. Hybridizations are carried out in the same solution and 5-20 X 
10 6 CpM 32 P- labeled probe is used. Filters are incubated in the hybridization mixture 
for 18-20 h at 55 °C, and then washed twice for 30 minutes at 60° C in a solution 
containing IX SSC and 0.1% SDS. Filters are blotted dry and exposed for 

20 autoradiography. Other conditions of moderate stringency which may be used are 
well-known in the art. Washing of filters is done at 37°C for 1 h in a solution 
containing 2X SSC, 0.1% SDS. 

For expression cloning (a technique commonly used in the art), an 
expression library is constructed. For example, mRNA is isolated from the cell type 

25 of interest, cDNA is made and ligated into an expression vector (e.g., a bacteriophage 
derivative) such that it is capable of being expressed by a host cell (e.g., a bacterium) 
into which it is then introduced. Various screening assays can then be used to select 
for the expressed gene product of interest based on the physical, chemical, or 
immunological properties of its expressed product. Such properties can be deduced 

30 from the properties of the corresponding orthologs from other species. 
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In another embodiment, polymerase chain reaction (PCR) can be used 
to amplify the desired sequence from a genomic or cDNA library. To isolate 
orthologous or paralogous genes from other species, one synthesizes several different 
degenerate primers, for use in PCR reactions. In a preferred aspect, the 
5 oligonucleotide primers represent at least part of the gene comprising known ortholog 
or paralog sequences of different species. It is also possible to vary the stringency of 
hybridization conditions used in priming the PCR reactions, to allow for greater or 
lesser degrees of nucleotide sequence similarity between the known nucleotide 
sequences and the nucleic acid homolog being isolated. 

10 Synthetic oligonucleotides may be utilized as primers to amplify by 

PCR sequences from a source (RNA or DNA), preferably a cDNA library, of potential 
interest. PCR can be carried out, e.g., by use of a Perkin-Elmer Cetus thermal cycler 
and a thermostable polymerase, e.g., Amplitaq (Perkin-Elmer), The nucleic acids 
being amplified can include mRNA or cDN A or genomic DNA from any eukaryotic 

1 5 species. After successful amplification of a segment of a the gene of interest, that 

segment may be molecularly cloned and sequenced, and utilized as a probe to isolate a 
complete cDNA or genomic clone. 

Once identified and isolated the gene of interest can then be inserted 
into an appropriate cloning vector for amplification and/or expression in a host. A 

* 

20 large number of vector-host systems known in the art may be used. Possible vectors 
include, but are not limited to, plasmids and modified viruses, but the vector system 
must be compatible with the host cell used. Such vectors include, but are not limited 
to, bacteriophages such as lambda derivatives, or plasmids such as pBR322 or pUC 
plasmid derivatives or the Bluescript vector (Stratagene). The insertion into a cloning 

25 vector can, for example, be accomplished by ligating the DNA fragment into a cloning 
vector which has complementary cohesive termini. 



6. EXAMPLE: USE OF SPECIALIZED DATABASES 

FOR IDENTIFICATION OF NOVEL GENES 
To test the method of using databases for gene discovery, protein 
30 sequence and domain/motif databases specific to two overlapping functional 
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groupings of proteins: (i) proteins known to be tumor suppressors, and (ii) proteins 
implicated in apoptosis in animals were developed. 

6.1 APOPTOSIS GENE DISCOV ERY METHOD 
5 Identification of a putative apoptosis-related human gene began with 

an identification of all genes in C elegans that contained either a POZ or kelch 
domain. A subset of these genes is shown in Figure 13. Hidden Markov Models 
(HMM) for the POZ and Kelch domains were built as follows. Starting with POZ and 
kelch sequences from the Drosophilia kelch protein (gi 1 577275) homologs were 

10 identified in other protein sequences using the BLASTP program. The resulting 

sequences showing significant similarity (e- value less than 0.001) were aligned using 
CLUSTAL W program and the alignments were used to build Hidden Markov Models 
with HMMER-2 package (Krogh et al., 1995, :http://hmmer.wustl.edu/). A computer 
printout listing of HMM models of tumor suppressors appears as a Microfiche H to 

15 the present specification. (See, http://hmmenwustl.edu; Chapter 2, which is 

incorporated by reference herein in its entirety, for a detailed description of HMM 
models) 

The resulting models were used to search through a database collection 
of C. elegans protein sequences. The domain structures of proteins having either a 

20 POZ or kelch domain were identified using existing collections of protein domains 
(e.g., see http://blocks.fhcrc.org^locks/blocks release.html, http://coot.embl- 
heidelberg.de/SMART/, http://www.motif.genome.ad.jp/). One of the unannotated 
protein-coding genes of C. elegans (corresponding protein accession number 
gi| 1 132541, see Figure 1 1) appeared to include a POZ domain, death domain, kinase 

25 domain, and heat repeat. A death domain is characteristic for the apoptosis system 
and a kinase domain indicates that the protein is likely to participate in 
phosphorylation of other proteins. The presence of these particular domains suggests 
that this protein is serving as a regulatory protein. 

Using the protein sequence of gi 1 1 132541, the database of human EST 

30 sequences was searched and a number of partial human cDNA sequences representing 
potential human orthologs or paralogs of the C.elegans gi| 1 132541 were identified. 
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The two closest human sequences, AA481214 and W51957, are depicted in Figure 
14 A. To determine whether the identified human sequences were orthologs or 
paralogs to the gi| 1 132541 gene of C. elegans, a gene tree (Saito and Nei, 1997, 
Molecular Biol. Evol. 4:406-425) was computed. The gene tree was generated using 
5 homologous genes identified with a BLASTP search against NCBI non-redundant 
database, using the human EST AA481214 sequence as a query. The resulting tree 
indicates that the identified human EST AA481214 represents a true ortholog of the 
Celegans gene gi 1 1 1 32541 (Figure 14B). The nucleotide sequence of the death 
domain protein is shown in Figure 14C, as well as the deduced amino acid sequence 
10 presented in Figure 14D. 

6.1.2 APOPTOSIS GENE DISCOVERY METHOD 
As a first step in identifying a novel gene involved in apoptosis, a 
comprehensive set of articles describing the system of apoptosis/programmed cell 

1 5 death in different species was compiled using the keyword "apoptosis". By analyzing 
the articles, information on regulatory pathways characterizing this system in different 
species, i.e., C. elegans, mouse, fruit fly, chicken, and human, was extracted. The 
regulatory information was stored as a collection of schemes produced in PowerPoint 
(Microsoft). Figure 4 shows a set of keywords defining proteins involved in apoptosis 

20 pathways. The keywords were used to generate a specialized sequence database, 
referred to as Apoptosis3, utilizing the PsiRetriever program for extraction of 
proteins from the all-inclusive non-redundant GenBank database (NCBI). Using 
program PsiRetriever, sequences from the non-redundant (NCBI) database of protein 
sequences, were retrieved and stored as a FASTA file. The FASTA file was then 

25 converted into binary blast database using program FORMATDB from the BLAST 
suit of programs. 

Genomic and cDNA sequences located in the region of human 
chromosome 13q were compared with the Apoptosis3 database using BLASTALL 
program from BLAST program complex. This region of the human genome is 

30 associated with Chronic Lymphocytic Leukemia (CLL). The comparison revealed 
significant similarity between a CLL region open reading frame and the mouse RPT1 
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protein (sp|P 1 5533 |RPT1) (Figure 13). Analysis of regulatory functions of RPT1 in 
the mouse reveals that this gene functions as a repressor of the interleukin 2 receptor 
(IL-2R) gene. When the RPT1 gene is knocked out, the regulatory effect is 
manifested as a block of the apoptotic pathway in T lymphocytes resulting in an 
5 accumulation of T lymphocytes in blood. This result is consistent with aberrations 
observed in CLL, namely abnormal accumulation of B-cells in the blood ( Trentin L. 
et al., 1997, Leuk. Lymphoma 27:35-42) and mutations in the human RPT1 gene play 
a role in development of CLL. 

6. 1 .3 EXAMPLE: A DISCOVERY OF A HUMAN ORTHOLOG OF THE 

1 0 MURINE MAX-INTERACTING TRANSCRIPTIONAL 

REPRESSOR : 

The family of Myc proto-oncogenes encodes a set of transcription 
factors implicated in regulation of cell proliferation, differentiation, transformation 
and apoptosis. C-Myc null mutations result in retarded growth and development of 

15 mouse embryos and are lethal by 9-10 day of gestation. In contrast, overexpression of 
Myc genes inhibits cell differentiation and leads to neoplastic transformation. 
Moreover, deregulation of Myc expression by retroviral transduction, chromosomal 
translocation or gene amplification is linked to a broad range of naturally occurring 
tumors in humans and other species. 

20 Another protein, called Max, is an obligatory heterodimeric partner for 

Myc proteins in mediating their function as activators of transcription during cell cycle 
progression, neoplastic transformation and programmed cell death (apoptosis). In 
order to make an active transcription factor the Myc proteins must form heterodimers 
with Max protein. This interaction with Max protein is necessary for specific binding 

25 of Myc with CACGTG box (or related E-boxes) on DNA and for activation of 
promoters located proximal to the binding sites. 

Besides the Myc family of transcription factors, the Max protein forms 
complexes with another family of so-called MAD proteins: Mxil, MAD1, MAD3 and 
MAD4. Whereas Myc:Max complexes activate transcription, MAD:Max complexes 

30 work in an opposite way repressing the transcription through the same E-box binding 
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sites and apparently antagonize Myc-mediated activation of the same set of target 
genes. 

During tissue development a shift from Myc\Max to MAD:Max 
complexes occurs coincidentally with the switch from cell proliferation to 
5 differentiation. The switch in heterocomplexes is thought to reflect a switch from 

activation to repression of common genes leading to cessation of proliferation, exiting 
the cell cycle and the beginning of cell differentiation. In differentiating neurons, 
primary keratinocytes, myeloid cell lines and probably other tissues the expression of 
different MAD'Max complexes appear in sequential order during the transition from 

10 cell proliferation to differentiation. The MAD3 expression appears first and it is 

restricted to proliferating cells prior to differentiation where it is co-expressed with 
two different member of Myc family, c-Myc or N-Afyc. Mxil transcripts are detected 
in proliferating and differentiating cells whereas MAD1 and MAD4 were confined to 
post-mitotic cells. Because Myc expression is not always downregulated in 

15 post-mitotic cells, co-expression of Myc and MAD genes may result in competition for 
Max heterodimers thus providing promoting or inhibitory effect on cell proliferation. 

The gene expression patterns, along with ability of Mad proteins to 
suppress A/yodependent transformation, are consistent with a potential function of 
Mad genes as tumor suppressors. This view is supported by the fact that allelic loss 

20 and mutations were detected at the Mxil locus in prostate cancers (Eagle et al., 1995 
Nat Genet 9:249-55). Cloning of the murine proteins Mad3 and Mad4 as well as their 
relation to Max signaling network was described by Hurlin (Hurlin PJ, et al., 1995, 
EMBO J. 14:5646-59) and Queva (Queva et al. 1998 Oncogene 16:967-977). Human 
orthologs of Mad4, Madl and Mxil are known. 

25 In this example, the discovery of an unknown human ortholog of Mad3 

protein found "in silico" by means of phylogenetic analysis of known mouse and 
human members of the Mad gene family and database searches is described. Since the 
function of murine Mad3 as a Mzx-interacting transcriptional repressor of 
jWyc-induced neoplastic transformation is well described, we can assign the same 

30 function to its human ortholog. 
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The gene tree shown in the Figure 20 was constructed in the following way. The 
protein sequences of known members of Mad gene family were extracted from 
GenBank database using NCBI Entrez keyword searches. The extracted sequences 
were aligned using multiple alignment program Clustalw running on Sun SPARC 
5 station. The quality of the multiple alignment was checked using program HitViewer 
Iterate (A. Rzhetsky, available upon request) and the redundant, non-homologous 
sequences as well as distant homologs from 5. cerevisiae, C elegans, D. 
melanogaster etc. were removed from the alignment. The refined set of sequences was 
realigned with Clustalw and a gene tree as presented in Figure 1 5 A was computed 

1 0 from the alignment using program NJBOOT (http://genome6.cpmc.columbia.edu // 
andrey) running on Sun SPARC station and viewed with program TreeView 
(http://genome6.cpmc.columbia.edu // andrey). 

The tree presented in Fig.l9A clearly shows the relationships between 
three known mouse genes and their two human homologs. Attempts to find a missing 

15 human ortholog of the mouse Mad3 gene in protein non-redundant database at NCBI 
using BLAST search did not identify any human homologs other than sequences that 
were already present on the tree, confirming the absence of a known human ortholog 
for Mad3 protein in the database. 

In order to identify a human ortholog of the Mad3 protein, a human 

20 dbEST at NCBI was searched with program TBLASTN using Mad3 protein sequence 
as a query. Two EST were identified and are shown in Figure 17A. 

Due to the nature of dbEST database this search produced only partial 
sequences of potential candidate genes. To obtain complete coding sequences 
(complete cds) of the genes, a search was conducted to obtain overlapping sequences 

25 in dbEST. The search for overlapping sequences was performed using the program 

Iterate with EST zs77e55.rl (gb|AA278224) serving as a query. The search returned a 
single overlapping sequence, namely HUMGSOO 12279 (dbj|C02407), thus indicating 
that the two EST sequences found during the initial TBLASTN search belong to the 
same gene. 

30 The complete sequence of the gene was assembled from the two ESTs 

using commercially available sequence assembly program SeqManll (DNASTAR Inc., 
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WI). The nucleotide sequence of the human Madi gene is presented in Figure 1 7B . The 
deduced amino acid sequence of the gene is presented in Figure 17C. The translated 
sequence consists of 206 amino acid residues 81% of which are identical to mouse Mad3 
protein. The alignment of human and mouse Mad3 proteins shown below was made using 
5 BLAST server at NCBI and is presented in Figure 1 7C. 

Multiple alignment of the new sequence with sequences of known Mad 
proteins was made using Clustalw and viewed with the HitViewer. A gene tree was 
computed from this alignment using NJBOOT. Multiple alignment of the new 
sequence with sequences of known Mad proteins (Figure 17C) along with its position 
10 on gene tree (Figure 18B) shows that this new human gene found by the approach 
described above belongs to the family of Mad proteins and is the ortholog of mouse 
Mad3. 

The present invention is not to be limited in scope by the specific 
embodiments described herein, which are intended as single illustrations of individual 

15 aspects of the invention, and functionally equivalent methods and components are 
within the scope of the invention. Indeed, various modifications of the invention, in 
addition to those shown and described herein will become apparent to those skilled in 
the art from the foregoing description and accompanying drawings. Such 
modifications are intended to fall within the scope of the appended claims. 

20 Various publications are cited herein, the contents of which are hereby 

incorporated by reference in their entireties. 
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CLAIMS 

1 . A method for identifying a novel nucleic acid molecule encoding 

a protein of interest comprising: 

(i) selecting a specific protein from a first species involved 
5 in a regulatory network of interest; 

(ii) identifying known proteins that act upstream and 
downstream in the regulatory network of interest with 
respect to the specific protein selected; 

(iii) constructing the regulatory network of interest from the 
10 proteins identified in step (ii); 

(iv) for each identified protein, select a domain or motif and 
search by homology for related proteins in a second 
species, wherein a related protein is defined as a protein 
having a homologous domain or motif; 

1 5 (v) producing a regulatory network for the second species, 

wherein said regulatory network incorporates the 

identified related proteins; 
(vi) comparing the regulatory network from the first species to 

the regulatory network of said second species; 
20 (v) identifying a protein present in a regulatory network for 

one species but absent in the regulatory network of the 

other species; and 
(vi) isolating a nucleic acid molecule encoding the protein 

identified in step (v) in the species in which it is missing. 

25 2. The method of Claim 1 wherein the nucleic acid molecule encodes 

human protein. 



3. The method of claim 1 wherein the related proteins are orthologs. 
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4. The method of claim 1 wherein the regulatory pathway is involved 

in apoptosis. 

5. The method of claim 1 wherein the specific protein from the first 

species is involved in tumor suppression. 

6. A method for identifying the affect of a gene knockout on a 
regulatory pathway comprising the following steps: 

(i) identification of the shortest non-oriented pathway 
connecting two gene products; 

(ii) assigning an initial sign value of "-" to the knockout since 
the knockout gene product is inactive; 

(iii) moving along the shortest pathway between the two gene 
products multiplying the sign with the sign of the next 
gene product in the pathway, wherein "-" stands for 
inhibition, "+" stands for induction or activation, and "0" 
stands for the lack of interaction between two proteins in 
the specified direction; and 

(iv) determining the final sign at the end of the pathway, 
wherein indicates inhibition and indicates 
induction or activation of the pathway. 

7. A method for identifying a novel nucleic acid molecule encoding 
a protein of interest comprising: 

(i) selecting a gene of interest and searching a database for 
homologous sequences; 

(ii) aligning the homologous sequences identified in step (j); 

(iii) constructing a gene tree using the sequence alignment; 

(iv) constructing a species tree; 
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(v) imputing the species tree and gene tree into an algorithm 
which integrates the species tree and the gene tree into a 
reconciled tree; and 

(vi) identifying orthologous genes present in one species but 

5 missing in another. 

8. The method of claim 7 wherein the following algorithm is used 
to integrate the species tree and the gene tree into a reconciled tree: 

(i) computing the similarity o(S gi JS 5 J) for each pair of 
interior nodes from trees T g and T s , 
1 0 (ii) finding the maximum o(S gi jS SJ ) ; 

(iii) saving S gi as a new cluster of orthologs, save {S gi } - {S sj } 
as a set of species that are likely to have gene of this kind 
(or lost it in evolution); 

(iv) eliminating S gj from T g ; T g : = T g \S g{ ; 

1 5 (v) repeating step (ii)-(iv) until T g is non-empty. 

9. A method for identifying a novel gene comprising the following 

steps: 

(i) defining a motif or domain composition of a gene of 
interest; 

20 (ii) searching for sequences which correspond to nucleotide 

sequences in an expression sequence tag database or other 
cDNA databases using a program such as BLAST and 
retrieving the identified sequences; 

(iii) searching additional databases for expressed sequence 
25 tags containing the domains and motifs characteristic for 

the gene of interest with Hidden Markov Model of 
domains and motifs identified in step (i); 

(iv) identifying nucleotide sequences comprising the gene of 
interest. 



WO 00/63687 



PCT/US00/10302 



53 



10. 



The method of claim 9 further comprising using each identified 
expression sequence tag to search sequence databases for 
overlapping sequences for the purpose of assembling longer 
overlapping stretches of DNA. 



5 



11. 



A method for extracting information on interactions between 



10 



biological entities 
(0 



(ii) 



from natural-language text data, comprising: 
parsing the text data to determine the grammatical structure of 
the text data ;and 

regularizing the parsed text data to form structured word terms. 



12. The method according to claim 1 1 , further comprising 
preprocessing the data prior to parsing, with preprocessing comprising the step of 
identifying biological entities. 

1 3 . The method according to claim 1 1 , further comprising referring 
15 to an additional parameter which is indicative of the degree to which subphrase 

parsing is to be carried out. 



further comprises segmenting the text data by sentences. 

1 5 . The method according to claim 1 1 , wherein said parsing step 
20 further comprises: 



14. The method according to claim 1 1 , wherein said parsing step 



segmenting the text data by sentences; and 

segmenting each of the sentences at identified words or phrases. 



1 6. The method according to claim 1 1 , wherein said parsing step 



further comprises: 



25 



segmenting the text data by sentences; and 
segmenting each of the sentences at a prefix. 
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17. The method according to claim 1 1 , wherein said parsing step 
further comprises skipping undefined words. 

18. The method according to claim 1 1 , wherein said parsing step 
further comprises: 

5 identifying one or more binary actions and their relationships; and 

identifying one or more arguments associated with the actions. 

1 9. The method according to claim 1 1 , further comprising 
performing error recovery when parsing of the text data is unsuccessful. 

20. The method according to claim 19, wherein said error recovery 
10 step comprises: 

segmenting the text data; and 

analyzing the segmented text data to achieve at least a partial parsing of 
the unsuccessfully parsed text data. 

2 1 . The method according to claim 1 1 , wherein said tagging step 

1 5 comprises providing the structured data component in a Standard Generalized Markup 
Language (SGML) compatible format. 

22. A computer system for extracting information on biological 
entities from natural-language text data, comprising: 

(i) means for parsing the natural-language text data; and 
20 (ii) means for regularizing the parsed text data to form structured 

word terms. 



25 



23. The system according to claim 22, further comprising means for 
preprocessing the data prior to parsing, with the preprocessing means comprising 
identifying biological entities. 
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24. The system according to claim 22, further comprising means for 
referring to an additional parameter which is indicative of the degree to which 
subphrase parsing is to be carried out. 

25. The system according to claim 22, wherein said parsing means 
5 further comprises means for segmenting the text data by sentences, 

26. The system according to claim 22, wherein said parsing means 
further comprises: 



27. The system according to claim 22, wherein said parsing means 
further comprises: 



means 



means 



for segmenting the text data by sentences; and 

for segmenting each of the sentences at identified words or 



10 phrases. 



means 



means 



for segmenting the text data by sentences; and 
for segmenting each of the sentences at a prefix. 



28. The system according to claim 22, wherein said parsing means 
further comprises means for skipping undefined words. 



29. The system according to claim 22, wherein said parsing means 



further comprises: 



means for identifying one or more binary actions and their relationships; 



20 and 



means for identifying one or more arguments associated with the actions. 



30. The system according to claim 22, further comprising means for 
performing error recovery when parsing of the text data is unsuccessful. 
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3 1 . The system according to claim 22, wherein said error recovery 
means comprises: 

means for segmenting the text data; and 

means for analyzing the segmented text data to achieve at least a partial 
5 parsing of the unsuccessfully parsed text data. 



32. The system according to claim 22, wherein said tagging means 
comprises means for providing the structured data component in a Standard 
Generalized Markup Language (SGML) compatible format. 
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p38 / p42 / ERK1 / p44 / ERK2 / SAFE / JNK / MEK / C-JTUN / MEF2D / ATFZ /emictamrto / ELK-1 / protein phosptou* 2A / 
rmf-1 / IH bete / TNF / FTK / Apof / p35 / ETS / C-Myc / XL* I receptor / NFOoipp* B / TNFR-1 / TRAIL / Apo-2L / 
DR4 / dcth rrocpinr / DR3 / DR2 / DR5 / DR1 / Ud / BMPR / BMP-* / TCF/ grim / bid / FAN / perforin / Fvl / F«s / DcRl 
/ decor receptor / w*»-l / NGF receptor / growth lector /RAR 
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>gi I 221 07 66 j gb |AA4 81214 t AA4 8121 4 aa34e02.rl NCI CGAP GCB1 Homo sapiens CDNA clone 
IMAGE: 815162 5* similar to WP:W07A12.4 CE03795 7, mRNA sequence (Homo sapiens] 
CATGGCTTCCTGGACACCAACCCTGCCATCCGGGAGCAGACGGTCAAGTCCATGCTGCTCCTGGCCCCAA 
AGCTGAACGAGGCCAACCTCAATGTGGAGCTGATGAAGCACTTTGCAC^ 

GGGCC CCAT CCG CT GCAACACCACAG TCT GC CTGG GCAAAAT C GGCTCCT ACCTCAGTGCTAGCACCAGA 
CACAG<K3TCCTTACCTCTGCCTTCAGCCGAGCCACTAGGGACCCGTTTGCACCGTCCCGGGTTGCGGGTG 
TCCT G<K»CT TT GCTGCCACCC^CAACCTCT ACTCAATG AACGACTGTG^ 
CTTGCGGTCTCACTGTAGATCCTGAGAAATCCGTGCGAGACC^VGGCCTTCAAGGCA 



>gi 1 134 9211 1 gb|W5l 957 JW519S7 zc45f01.rl Soares_seneseent_f ibroblssts_NbHSF Homo 

sapiens cDNA clone IMAGE: 325273 5 f , mRNA sequence [Homo sapiens] 

CCTTCGAGTTCGGCAATGCTGGGGCCGTTGTCCTCACGCCCCTCTTCAAGGTGGGCAAGTTCCTGAGCGC 

T GAGG AGT AT CAGCAGAAGATCATCCCTGTGGTGGT CAAG ATGTTCTCAT CCACTGACCGGGCCATGCGC 

ATCCGN CTCCT GCAGC AG AT GGAGCAGTTCATCCAGT ACCT TGACGAGCCAACAGTCAACACCCAGAT CT 

TCCCCCACGTCGTACATGGCTTCCTGGACACCAACCCTGCCATCCGGGAGCAGACGGTCAAGTCC^TGCT 

GCTCCT GGCCCCAAAGCTGAACGAGGCCAACCTC AATGT GG AGCT GATGAAGCACTTTGCACGGCTACAG 

GCCAAGGATGAACAGGGCCCC^TCCGCTGCAACACCAC^ 

GTGCTAGCACCAGACACAGGGTCCTTACCTCTG 
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S.cervisiae_S60992 

Nicotiana_tabacum_e244568 



0.15 



WO 00/63687 



16/23 



PCT/US00/10302 



BASE COVNT 
ORIGIN 



405 a 



545. c 



493 g 



278 t 



6 others 



1 cagccgaagc aiugcaaaaat tcttccagga gctgagcaag . agcctggacg Catfcccctga 
61 ggayttctgt cggcacaagg tgctgcccca gctgctgacc gccttcgagt tcggcaatgc 
121 tggggccgtt gtcctcacgc ccctcttcaa ggtgggcaag ttcctgagcg ctgaggagta 
181 tcagcagaag atcatccctg tggtggtcaa gatgttctca tccactgacc gggccatgcg 
241 catccgccrtc ctgcagcaga tggagcagtt catccagtac cttgacgagc- caacagtcaa 
301 cacccagatc ttcccccacg tcgtacatgg cttcctggac accaaccctg ccatccggga 
361 gcagacggtc aagtccatgc tgctcctggc cccaaagctg aacgaggcca acctcaatgt 
421 ggagctgatg aagcactttg cacggctaca ggccaaggat gaacagggcc ccatccgctg 
4 81 caacaccaca gtctgcctgg gcaaaatcgg ctccfcacctc agtgctagca ccagacacag 
541 ggtccttacc tctgccttca gccgagccac tagggacccg tttgcaccgt cccgggttgc 
601 gggtgtcctg ggctttgctg ccacccacaa cctctactca atgaacgact gtgcccagaa 
661 gatcctgcct gtgctctgcg gtctcactgt agatcctgag aaatccgtgc gagaccaggc 
721 cttcaaggcm wttcggagct tcctgtccaa attggagtct gtgtcggagg acccgaccca 
781 gctggaggaa gtggagaagg .atgtccatgc agcctccagc cctggcatgg gaggagccgc 
841 agctagctgg gcaggctggg cgtgaccggg gtctcctcac tcacctccaa gctgatccgt 
901 tcgcacccaa ccactgcccc aacagaaacc aacattcccc aaagacccac gcctgaagga 
961 gttcctgccc cagcccccac ccctgttcct gccaccccta caacctcagg ccactgggag 
1021 acgcaggagg aggacaagga cacagcagag gacagcagca ctgctgacag atgggacgac 
1081 gaagactggg gcagcctgga gcaggaggcc gagtctgtgc tggcccagca ggacgactgg 
1141 agcaccgggg gccaagtgag ccgtgctagt caggtcagca actcbgacca caaatcctcc 
1201 aaatccccag agtccgactg gagcagctgg gaarctgagg gctcctggga acagggcrtgg 
1261 caggagccaa gctcccagga gccacctyct : gacggtacac ggctggccag cgagtataac 
1321 tggggtggcc cagagtccag cgacaagggc gaccccttcg ctaccctcftc tgcacgtccc 
1381 agcacccagc cgaggccaga ctcttggggt gaggacaact gggagggcct cgagactgac 
1441 agtcgacagg tcaaggctga gctggcccgg aagaagcgcg aggagcggcg gcgggagatg 
1501 gaggccaaac gcgccgagag gaaggtgcca. agggccccat gaagctggga gcccggaagc 
1561 tggactgaac cgtggcggtg gcccttcccg gctgcggaga <?cccgcccca cagatgtatt 
1621 tattgtacaa accatgtgag cccggccgcc cagccaggcc atctcacgtg tacataatca 
1681 gagccacaat aaatrtctatt tcacaaaaaa aaaaaaaaaa aaaaaaa 



// 
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>sp|P15533|RPTl JMOUSE DOWN REGULATORY PROTEIN 
OF INTERLEUKIN 2 RECEPTOR (J03776) rpt-lr [Mus 
musculus] Length « 353 

Score « 92.0 bits (237), Expect = 6e-20 

Query 194 VMELLEEDLTCPICCSLFDDPRVLPC5HNFCKKCLEGII^GSVIWSMWRPAPFKCI>TCRK 373 

V+E+++E++TCPIC 1» +P C+H+PC+ C+ E S RN+ CP CR 

Sbjct 5 VLEMIKEEVTCPICLEIAKEPVSAIXZNHSrCRACITLNYE-SNRNT DGKGNCPVCRV 60 

Query 374 ETSATGIN5LQVNYSLKGIVEKYNKIKI5F KMPVCKGHMGQPLNI FCLTDMQLICG 541 

- +L+ N + IVE+ K P K+ +C H G+ L +FC DM +IC 
Sbjct 61 PYP FGNLRPNLHVANrVERLKGFKS I PEEEQKVNX CAQH-GEKLRLFCRKDMMVICW 116 

Query 542 I C ATRGEHTKHVFC S I EDAYAQERDAFESLFQS F ETWRRGDALSR1DTMETSK 700 

+C EH H IE+ + ++ + "W+ L R+D 
Sbjct 117 LCERSQEHRGHQTALI EEVDQE YKEKI^GALWKLMKXAICI CDEWQDDLQLQRVDW 171 

Query 701 RKSLQLMTKDSDKVKEFFEKLQHTLDQKKNEILSDFETMKLAV^ 862 

+Q+ + + V+ F+ L+ •LD K+NE L + K VM+ + N+L 
Sbjct 172 ENQIQI NVENVQRQFKGIJU3L1J3S KZNEEIX}KLKKEKKEVKEKLE ES ENEL. 222 



Homology covers ring finger, B-box and the beginning of coiled coil domain 
in the CLL ring finger protein 



Figure '5 
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Activated CD4* T-cells 

Rptl (represses expression of IL-2 receptor) 



IL-2 receptor 



> normal expression of Bcl2 



IL-2, EL-15 



normal apoptosis 



When rptl is knocked out: 

IL-2, IL- 15 

I 

IL-2 receptol ► overexpression of Bcl2 

apoptosis 
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TBIASTH 2.0.8 (Jvi»O5-l«90] 

Reference : 

Altschul, Stephen F., Thomas L. Madden, Alejandro A. Schaffer, 
Jinghui Zhang, Zheng Zhang, Webb Miller, and David J. Lipman (199*7), 
"Gapped BLAST and PSI -BLAST: a new generation of protein database search 
programs". Nucleic Acids Res. 25:3389-3402* 

Query- gi I 2137498|Mad3m 
(205 letters) 

gblAA27B224 1AA278224 *$77eQ5.rl NCI CGAP_GCB1 Homo sapiens cDNA clone IMAGE:703520 5* 
similar to TR:G1184157 gI184157 MAX- INTERACTING 
TRANSCRIPTIONAL REPRESSOR. ; 
Length * 430 

Score - 209 bits (526) , Expect • le-53 

Identities - 104/124 (83%), Positives - 116/124 (92%), Gaps - 1/124 (0*) 
frame -42 

Query: 1 ME PVASNI QVLLQAAE FLERREREAEHGY ASLC PHHS PGTVCRRRKP PLQAPGALN SGRS 60 

MEP+ASNIQVXLQAAEFLERREREAEHGYASLCPH SPG + RR+K P QAPGA +SGRS 
SbjCt: 56 ME PLASNI QVLLQAAE FLER RE REAEHG Y ASLC PHRS PG PI HRRKKR P PQA PGAQD SGRS 235 

Query: 61 VHNELEKRRRAQLKRCLEQLROOMPLGVDCTRYTTLSLL-RARVHIQKLEEO^^ 119 

VHNELEKRRRAQLKRCLE+ L+QQMPLG DC RYTTLSLL RAR+HI QKLE+QEO+AR+LK 
SbjCt: 236 WWELEKRRRAQLKRCLERLKQC^PLGGDC^YTTLSLLRRARMHIQKLEDQEQRARQLK 415 

Query: 120 EKLRS 124 
E+LR4 

SbjCt: 416 ERLRT 430 

dbj 1C02407(C02407 HUMGS0012279, Human Gene Signature, 3'-directed cDNA sequence. 
Length - 34 8 

Score « 97.5 bits (239), Expect - 6e-20 
Identities - 51/63 (80%), Positives - 56/63 (87%) 
Frame ■ +3 

Query: 125 KQQS LQQQLEQLQGL PGARE RERLRADSL DSSGLSSERSDS DQE DLEVDVENLV FGTETE 1B4 

KQQSLQ+ QL+GL GA ERERLRADSLDSSGLSSERSDSDQE+LEVDVE+LVETG E E 
SbJct: 4 5 KQOSLQRXWMQLRGIAGAAERERLRADSLDSSGLSSERSDSDQEELEVDVESLVFGGEAE 224 

Query: 185 LLQ 187 
LL+ 

Sbjct: 225 LLR 233 
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BASE COUNT 130 a 234 c 258 g 106 t 5 others 

ORIGIN 

1 cagccgcttg ctccggccgg caccctaggc cgcagtccgc caggctgtcg ccgacatgga 
61 acccttggcc agcaacatcc aggtcctgct gcaggcggcc gagttcctgg agcgccgtga 
121 gagagaggcc gagcatggtt atgcgtccct gtgcccgcat cgcagtccag gccccatcca 
181 caggaggaag aagcgacccc cccaggctcc tggcgcgcag gacagcgggc ggtcagtgca 
241 caatgaactg gagaagcgca ggagggccca gttgaagcgg tgcctggagc ggcfcgaagca 
301 gcagatgccc ctgggcggcg actgtgcccg gta caeca eg ctgagcctgc tgcgccgtgc 
361 caggatgeae atccagaagc tggaggatca ggagcagegg gcccgacagc tcaaggagag 
421 gctgcgcaca aagcagcaga gcctgcagcg gcantggatg cagctccggg ggfctggcagg 
481 ngeggecgag egggagegne tgcgggcgga eagtctggac tcetcaggcc tctcctctga 
541 gcgctcagac tcagaccaag aggagctgga ggtggatgtg gagagectgg tgtttggggg 
601 tgaggecgag ctgctgcggg gcttcgtcgc eggecaggag cacagctacrt cgcacgtcgg 
661 cggcgcctgg ctatgatgtt cctcacccan ggcgggcctc tgccctctta cfccgttgccc 
721 aagcccactt tnc 
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>Afa43b(Putative) 

MEP1ASKIQVLLQAAEFLERREREAEHGYASLCPHRSPGPIHRRKKRPPQAPGAQDSGRSVHNELEKRRRAQLK 
RCLERLKQQMPLGGDC^YTTLSLLRRARMHIQKLEDQEQR 

LRADSLDSSGLSSERSDSDOEELEVDVESLVFGGEAELLRGFVAGOEHSYSHVGGAWL 



gil2S068881MAIta 
gil729978IR«>h 
gi 12792362 IM«d4h 
gil2137499IM»d4» 
gi I 21374 99 IHad3v 
H*d3h Putativa 



KATAVC*HIQLL:Jy^/l.ERAEREAXK;Y^ 
HUUVR«JIJ(*1LLEAADYLERW:R£A£rc^ 

HELHSLL3 LLEAAEYLERRDREAEHCYASVLPFlXroFAAEKTKAAGIA'RKAP^ 

HE1J<5LLLLLEAAEYLERRDR£>EHGYASMLPFDGDFARXJCTKTAG^ 

-HEPVASH I QVLLQAAE F LERRE RE AERGY AS LCPHHS PGTVCRRRK P P LQAPGALMS GRSVHNELEKRRRAQLXRCLEQLRQQMPLGVt<rTRYTTLS LL 
-HE P LAS K I CvT,LGAAI F RAE REAE BG7 AS LCP MRS PGP 1H RRXK RP PQAPG AQDS G R5 VHJfE LEXRRRAQUCRCIJE RLKQ04 P1GGDCAR YTTLS LL 



gi | 2S06BBB IMADn 
911729978 IKADh 
gi|2792362 IK*d*h 
gi 12137499 lK»dl* 
gil21374»8tH«43» 
H«d3h Put«tiv« 



TKAKLH1 KKU3)O)RKAVH0J DOWttQRHLKRALEKLGAERTR KOS VC-S WSSER50SDREEUJV^T)VX^TM>^TDYUIODl^ 

TKAKLHIKKUDCDRKAVHQIOQLOREORHLKRQLEKLCIERIR KDS 1G-STVSSERSDSDRE EIDVDVESTOYLTCDLDWSSSS 

! KKLEEQORRALS IKEOM)QEKRFLXRRLEQU« VQSVIR VRTD5TG-SAV5TD — PSEQE * VDl BCME FGPG E - 

K-AKHHI KKLEEOORRALS IKEOI.QREHRFLXRRLEQLS VQ5VU VRTD5TG-5AVSTD — DSECE VDl EGWEFGPOEUttAGS - 

R-AAVHlQKLEEWQQ*RRI4aKlil5KQQSl^LEQL^l.WSAR£P^R^ DLSVDVEHLVFG-TET8X.LQS F 

AAAAMHICKLEDOEORARGIJCERLRTKQOSWRXwto^ EUFVUVES LVFC-GEAELLRGF 



gj 1 25068 88 IMAEtt 
giT?29978IKADh 
gi 12792362 |Mad4h 
gi|2137499lKad4B 
gil 2137498 IttacUm 
H*d3h Putativt 



VSDSDERGSMQS LG-S DEGY5S ArvTCRAKLQCXSHKAG jjG L 
VS DS DERGSMQS 1G-S DEG YSS TStKR IKLQDS HKACLC L 

SSDADDHYS LOS GTGGOSG FG PHCRR1GRPAL5 

5SDADDKYSl>OSSGCSDSSYCaFCRRPGCPGLS— 

SAGRERSYSHSTCAtfL < 

VAGOEHSYSHNW^WL — . 
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% lexsemsub.pl 

% lexsemsub. pat 

% revised March 17, 2 000 

% LEXICON OF SUBSTANCES AND STRUCTURES 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
: -multifile (phrase/5) . 
: -multifile (wdef/3 ) . 
: -unknown (_, fail) . 

phrase ( 1 [ 1 , protein, ['[', gamma, ' ]', ' -', aminobutyric , acid, a], 1 G A 
BAA 1 , r) . % ? 

phrase ('[', smallmolecule , [ 1 [ ' , zeta, ' ] ' , 1 , subunit] , 1 [zeta]l subu 
nit ' , r) . % ? 

phrase (116 , protein, [116, 1 - ' , kd, fyn, ' - ' , associated, protein] , 1 116-k 
D Fyn-associated protein', r) . 

phrase(116, protein, [116 , 1 -', kd, protein] , '116-kd protein', r), 

phrase (3 ,protein, [3 , 1 - ' , kinase, 1 - 1 , akt] , ' 3-kinase-Akt ' , r) . 

phrase (ability, affirmation, [ability, to] , [] , r) . 

phrase (age, protein, [age, protein, kinases], 1 AGC ' , r) . 

phrase (akt , protein, [akt, mutant], 'Akt mutant 1 , r) . 

phrase (alternative, substance , [alternative, ntf] , 'alternative NTF 1 , r 

) • 

phrase (antibody, protein, [antibody, to, phosphotyrosine] , 'anti-phosp 
hotyrosine 1 , r) . ■• 

phrase (antigen, complex, [ant igen, receptor] , 'antigen receptor' ,r) . 
phrase (ap, protein, [ap, 1 - ' , 1] , ' AP-1 1 , r) . 

phrase (aspargine, site, [aspargine, ' - 1 , 141] , 1 aspargine-141 ' , r) . 
phrase (b, cell, [b,cell], 'B cell', r) . 
phrase (b, cell, [b, cells], 'B cell', r) . 

phrase(b, species, [b, lymphoblastoid, cells] , ! B lymphoblastoid cell 
s 1 , r) . 

phrase (b, cell , [b, lymphoblastoid, cells] , 'B lymphoblastoid cells 1 , r 
) . 

phrase (b7, protein, [b7 , 1 - ■ , r 1 1 ] , , B7-l , ,r). 
phrase (bcl , protein, [bcl, 1 - ' , 2] , ' Bcl-2 ' , r) , 
phrase (c, protein, [c, 1 - 1 , jun] , 'c-Jun' ,r) . 
phrase (camk, protein, [camk, iv] , ' CaMK IV, r). 
phrase (casp, protein, [casp, ' - ' , 3] , ' caspase-3 1 , r) . 

phrase (caspase, protein, [caspase, ' -', 3 , family , protease] , 'caspase-3 

family protease' ,r) , 
phrase (caspase, protein, [caspase, ' - 1 , 3 , precursor] , 1 caspase-3 precur 
sor 1 , r) . 

phrase (caspase, protein, [caspase , ' - ' , 3 ] , ' caspase-3 1 , r) . 
phrase (caspase, protein, [caspase , - , 3] , ' caspase-3 ' , r) . 
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phrase (caspase, protein, [caspase, » - ' , 6] , 1 caspase -6 1 , r) . 

phrase (caspase, protein, [caspase, ■ - 1 , 7] , ' caspase- 7 • , r) . 

phrase (catalytic, domain, [catalytic, domain], 'catalytic domain', 

r) . " 

phrase (cleavage, site, [cleavage, site] , 'cleavage site',r) . 

phrase (cleavage, substance, [cleavage , products] , 'cleavage products' , 

r) . 

phrase (cooh, substance, [cooh, ' - ' , terminal , fragment] , ' COOH- termina 
1 fragment ' , r) . 

phrase (crk, protein, [crk, proteins] , 'crk proteins ', rO . 

phrase (crkl, complex, [crkl ,'-', c3g , complex] ,' crkl -c3g complex ' , r) . 

phrase (dcp, protein, [dcp, - , 1] , 'DCP-1 r) . 

phrase (did, negation, [did, not], not, r) . 

phrase (ebv, species, 1 Epstein-Barr virus' ,r) . 

phrase (epstein, species, [epstein, ' - « , barr, virus] , ■ Epstein-Barr vi 
rus ' , r) . 

phrase (familial , disease, [familial , alzheimer, ' 1 1 ' ,s, disease] , ' f amil 
ial Al zheimer • ■ ' * s disease ' , r) . 

phrase(gene, gene, [gene, encoding, interleukin, '-',2], 'gene encodin 
g interleukin-2 1 , r) . 

phrase (gst, protein, [gst, 1 - » , 1 f yn 1 , ' - ' , sh2] , ' GST-Fyn- SH2 1 , r) . 
phrase (gst, protein, [gst , ' - ' , ' fyn' , ' - ' , sh3] , ' GST-Fyn-SH3 1 , r ) . 
phrase(gtp, complex, [cjtp , exchange, of , rapl] , 'GTP exchange of Rapl 1 , 

phrase (guanidine, protein, [guanidine, nucleotide, ' - 1 , releasing, f ac 
tor,c3g] , 'guanidine nucleotide-releasing factor C2G~ f ,r) . 
phrase (guanidine, smallmolecule, [guanidine, nucleotide] , 'guanidine 
nucleotide' ,r) . 

phrase (guanosine, smallmolecule, [guanosine, triphosphate] , 'guanosin 
e triphosphate 1 , r) . 

phrase (guanosine, smallmolecule, [guanosine , diphosphate] , 'guanosine 
diphosphate' f r) . 

phrase (h4, cell, [h4, cell, line] , 'H4 cell line',r). 

phrase (h4, cell, [h4 , human, neuroglioma, cells] , ' H4 , human, neuroglioma 
, cells ' , r) . 

phrase (ha, protein, [ha, '-', '.[«, delta, '] *,phpkb] , ' HA- [Delta] PHPK 
B' , r) . 

phrase (hla, protein, [hla, ' - 1 , dr7] , 1 HLA-DR7 ' , r) . 

phrased, protein, [i, '[«, kappa, ']', b, ■ -','[', beta, ']»] , ' I [ka 
ppa] B- [beta] 1 , r) . 

phrase (i, protein, [i, '[', kappa, 1 ]', b ,' -','[', alpha, ']'] , 1 I [kap 
pa] B- [alpha] 1 , r) . 

phrase (i , protein, [i, '[', kappa, ' ]',b], ' I [kappa] B ' , r) . 
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phrase (ice, protein, [ice, 1 / 1 , ced, t - ! , 3] , 1 ICE/Ced-3 1 , r) . 

phrase (il, gene, [il , 1 - ' , 2 , gene] , 'gene encoding interleukin-2 ' , r 

) - 

phrase (il , protein, [il , 1 - 1 , 2 ] , ' interleukin-2 1 , r) . 

phrase (in, interm, [in, the, case, of] , [] , r) . 

phrase (in, state, [in, the, anergic, state] , inactive, r) . 

phrase (inducible, cell, [inducible , h4 , cell] , 'inducible H4 cell f ,r 

) . 

phrase (interleukin, protein, ['interleukin, 1 - 1 , 2] , r) . 

phrase (interleukin, protein, [interleukin, ' 3], ' interleukin-3 

',r) . 

phrase (interleukin, protein, [interleukin, ' - 1 ,1, beta, converting, enzy 

me], 1 interleukin-1 beta converting enzyme 1 , r). 

phrase (jurkat, cell, [jurkat, cell], 'Jurkat cell', r) . 

phrase (jurkat , cell, [jurkat, cells], 'Jurkat cell', r) . . 

phrase (kif 3a, protein, [kif3a, '/' ,3,b] , 'KIF3A/3B' , r) , 

phrase (lbl, cell, [lbl , ' - • , drf , cells], 1 LBL-DR7 cells', r). 

phrase (lbl, cell, [lbl, ' - • ,dr7, cells] , 1 LBL-DR7 cells ' , r) . 

phrasedet, protein, [let , 1 - • , 23] , 'Let-23 1 , r) . 

phrase (may, probability, [may, be] , possible, r) . 

phrase(myc, protein, [myc, p70s6kd3e] , ' Myc-p70s6kD3E 1 , r) . 

phrase (myc, protein, [myc, 1 - 1 , pdkl] , 'Myc-PDKl 1 ,r> . 

phrase (myc, protein, [myc, ' - 1 ,p70s6k] , , Myc-p70s6k' , r) . 

phrase (myc , protein, [myc, ' - 1 , p70s6ke3 89d3e] , ' Myc-p70s6kE389D3E ' , r) 

phrase (myr , protein, [myr , ' - 1 , akt] , ' Myr-Akt 1 , r) . 

phrase (n, protein, [n, 1 - ' , methyl , ' - 1 , d, 1 - 1 , aspartate, receptor] , 1 N 
MDAR 1 , r) . 

phrase (n, protein, [n, ' - ' , methyl, ' - ' , d, ' - ' , aspartate] , ' NMDA 1 ) . 
phrase (native, cell, [native, h4 , cell] , ' native H4 cellar), 
phrase (nf , protein, [nf ,'-','[• # kappa, ■ ] ■ ,b] , »NF- [kappa] B 1 , r) . 
phrase (nh2,. site, [nh2 , 1 - 1 , terminal] , 1 NH2- terminal 1 , r) . 
phrase (nh2 , substance, [nh2, terminal , fragment] , 'NH2-terminal fr 
agment 1 , r) . 

phrase (nih, cell, [nih, ' - 1 , 3 , t3 , fibroblasts] , 'NIH-3T3 fibroblasts 1 
i r ) • 

phrase (nih, cell, [nih, 1 - ' , 1 3t3 ■ , fibroblasts] , 'NIH-3T3 fibroblasts 1 
,r) . 

phrase (normal , substance, [normal, ntf] , 'normal NTF ' , r) . 

phrase (nuclear, protein, [nuclear, factor, kappa, b] , ' NF- [kappa] B ' 

, r) . 

phrase (plBOGlued, protein, [plSOGlued, - , arpl] , ' pl50Glued-Arpl ' , r ) . 
phrase (phosphate, phosphorylate2 , [phosphate, incorporated, into] , 
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phosphorylate , r) . 

phrase (phosphatidyl inositol , smallmolecule, [phosphatidylinositol , 1 
, , , , / 4, , / , ,5, 1 -', triphosphate] , 'phosphatidylinositol 1,4, 5-tripho 
sphate 1 , r) . 

phrase (phosphoinosi tide, protein, [phosphoinositide, 1 - ' , dependent, 

protein, kinase], ' PDKl',r). 
phrase (phospholipase, protein, [phospholipase, c, '-',1] , 'phospholip 
ase C- 1 ' , r) . 

phrase (poly, protein, [poly, * { ' , adp, ■ - 1 , ribose, 1 ) 1 , polymerase] , 1 poly 
(ADP-ribose) polymerase' , r) . 

phrase (polyvinylidene, structure, [polyvinyl idene, difluoride, memb 
ranes] , 'polyvinylidene difluoride membranes' , r) , 
• phrase (presenilin, protein, [presenilin, 1] , 'presenilin 1 1 , r) . 
phrase (presenilin, protein, [presenilin, 2] , ' presenilin, 2 • , r) . 
phrase (productively, state, [productively , stimulated] , active, r) * 
phrase (protein, protein, [protein, tyrosine , kinase] , 'protein tyrosi 
ne kinase ' , r) . 

phrase (protein, protein, [protein, kinase, c] , 'protein kinase C',r) . 
phrase (ps2 , substance, [ps2, '-',ctf] , 'presenilin 2 COOH- terminal fra 
gment 1 , r) . 

* phrase (ps2 , substance, [ps2 , cleavage, fragment] , 'presenilin 2 cleava 
ge fragment 1 , r) . 

phrase (pvdf, structure, [pvdf, membranes] , 'polyvinylidene difluori 
de membranes ' , r) . 

phrase (raf, protein, [raf , ' - ' , 1] , 'Raf-l 1 , r) . 
phrase (raf , protein, [raf, ■ - ' , 1] , 'Raf-l ' , r) . 
phrase (rapl , complex, [rapl , ' - ' , gtp] , ' Rapl-GTP ' , r) . 
phrase (requirement , need2, [requirement, for], need,r). 
phrase (ser, smallmolecule, [ser, 19], ' Ser 19 ',r). 
phrase (ser, smallmolecule, [ser, 23], 'Ser 23',r). 

phrase (serine, substance, [serine, residues], 'serine residues', r 
) . 

phrase(src, domain, [src, homology, 2], 'Src homology 2 ' , r) . 
phrase(src, domain, [src, homology, 3], 'Src homology 3 1 , r) . 
phrase (srebp, protein, [srebp, ' - 1 , 1] , ' sterol -regulatory element bin 
ding protein 1 1 , r) . 

phrase (srebp, protein, [srebp, '-', 2] , ' sterol -regulatory element bin 
ding protein 2' ,r) . 

phrase (sterol , protein, [sterol, ' - ' , regulatory , element , binding, prote 
in,l] , 1 sterol-regulatory element binding protein 1' ,r) . 
phrase (sterol, protein, [sterol, *-' , regulatory , element , binding, prote 
in, 2] , 'sterol-regulatory element binding protein 2' ,r) . 
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phrase (t, cell, [t, ■ - • , dr7] , 1 t-DR7 1 , r) . 

phrase (t, cell, [t, ' -' ,drt, ■ / 1 , b7 , ' - ' , 1] , » t-DR7/B7-l 1 ,r) . 
phrase(t, cell, [t,cell], ! T cell ! ,r). 
phrase(t, cell, [t, cells], 'T cellar). 

phrase (t , complex, [t , ' - ' , cell , receptor] , ' T-cell receptor 1 , r) . 

phrase (t , cell, [t , ' - * , dr7 , cells] , ' t-DR7 cells ' , r) . 

phrase (t, cell, [t , » - » , dr7, 1 / ' , b7 , ' - 1 , 1] , ■ t-DR7/B7-l ■ , r) . 

phrase (t , complex, [t, cell , antigen, receptor] , { T-cell antigen red 

eptor 1 , r) . 

phrase (threonine, aminoacid, [threonine, 229], 'threonine 229', r) 



phrase (transcription, ' protein, [transcription, factor] , 'transcript 
ion factor ' , r) . 

phrase (trypan, smallmolecule, 'trypan blue',r) . 

phrase (wt, protein, [wt, akt] , 'WT Akt',r) . 

phrase ( zap, protein, [zap, * - 1 , 70] , 1 ZAP-70 ' , r) . 

phrase (zdevd, smallmolecule, [zdevd, ' - 1 , fmk] , 1 zDEVD-f mk ' , r) . 

phrase (il , protein, [il,'-'^],' interleukin-3 ' , r) . 



wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 



ab, complex, antibody) . 

actin, protein, actin) . 

activated, state, active) . 

active, state, active) . 

ad, disease, 'Alzheimer' ■ 1 ' s disease ' ) . 

age, protein, 1 AGC ' ) . 

akt, protein, 'AKT') . 

anergic, state, inactive) . 

anergic , state, inactive) . 

anergy, state, inactive) . 

antibody , complex, antibody) . 

antigen, substance, antigen) . 

aop, protein, 'Aop'). 

apoptosis, process, apoptosis) . 

bad, protein, 'BAD'). 

c3g, protein, 1 C3G 1 ) . 

' ca2 + 1 , smallmolecule, ' Ca2+ ' ) . 

cas, protein, 'Cas'). 

caspase, protein, caspase) . 

caspase, protein, caspase) . 

cbl, protein, 'Cbl'). 

ccrsrh, protein, ' CCRSrh') . 

cd28, protein, f CD28'). 

cells, structure, cell) . 

cholesterol, smallmolecule, cholesterol) 
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wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 
wdef 



cpp32, protein, 'CPP32') . 
crkl, protein, ' CrkL'). 

ctf , substance, 1 COOH- terminal fragment' ) . 

cytokine, smallmolecule, cytokine) . 

cytosol, structure, cytosol) . 

djnk,protein, 1 D JNK 1 ) . 

djun, protein, 'DJun' ) . 

dynamitin, protein, dynamitin) . 

erk, protein, ' ERK 1 ). 

eto, smallmolecule, 1 ETO ' ) . 

etoposide, smallmolecule, etoposide) . 

fad, disease, 1 familial Alzheimer ' 1 '' s disease') 

fyn, protein, f Fyn'). 

gdp, smallmolecule, 'GDP') . 

gelsolin, protein, gelsolin) . 

gpl20, protein, 'gpl20') . 

grb2, protein, 'Grb2'). 

gst, protein, 'glutathione S- transferase 1 ) . 

gtp, smallmolecule, ' GTP 1 ) . 

hsp70 , protein, 1 HSP70 1 ) . 

human, species, human) . 

ikk, protein, 'IKK'). 

inactivated, state, inactive) . 

inactive, state, inactive) . 

jnk, protein, 1 JNK' ) . 

jnk, protein, 'JNK'). 

jnk2, protein, 1 JNK2 1 ) . 

kap3 , protein, kap3 ) . 

kdakt , protein, 1 KDAkt 1 ) . 

kinase, protein, kinase) . 

kinectin, protein, kinectin) . 

klc, protein, klc) . 

lamin, protein, lamin) . 

myosins, protein, myosins) . 

nmdar, protein, 1 NMDAR 1 ) * 

nmdar2b, protein, 1 NMDAR2B 1 ) . 

ntf, substance, 1 NH2 -terminal fragment') . 

p70s6k, protein, p70s6k) . 

p78s6k, protein, p78s6k) . 

parp, protein, 1 poly (ADP-ribose) polymerase ! ) . 
pdkl, protein, ' PDKl 1 ) . 
peptides, protein, peptide) . 
pkb, protein, 1 PKB 1 ) . 
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